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I. INTRODUCTION 

A. The physics context 

The theoretical precision of observables like the run- 
ning of the Quantum Electrodynamic (QED) fine struc- 
ture constant a(s) or the anomalous magnetic moment 
of the muon is limited by second-order loop effects from 
hadronic vacuum polarization (VP). Theoretical calcula- 
tions are related to hadronic production rates in e + e~ 
annihilation via dispersion relations. As perturbative 
Quantum Chromodynamic theory fails in the energy 
regions where resonances occur, measurements of the 
e+e~ — > hadrons cross section are necessary to evaluate 
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the dispersion integrals. Of particular interest is the con- 
tribution a^ ad to the muon magnetic moment anomaly 
a M , which requires data in a region dominated by the 
process e + e~ — > 7r + 7r~(7). The accuracy of the theoret- 
ical prediction for is linked to the advances in e + e~ 
measurements. A discrepancy of roughly 3 standard de- 
viations (<t) including systematic uncertainties between 
the measured [l[ and predicted 0-0] values of a M per- 
sisted for years before the results of this analysis became 
available [7|, possibly hinting at new physics. An inde- 
pendent approach using r decay data leads to a smaller 
difference of 1.8cr @ in the same direction, with enlarged 
systematic uncertainties due to isospin-breaking correc- 
tions. 

The kernel in the integrals involved in vacuum polar- 
ization calculations strongly emphasizes the low-energy 
part of the spectrum. 73% of the lowest-order hadronic 
contribution is provided by the 7r + 7r _ (7) final state, 
and about 60% of its total uncertainty stems from that 
mode Q. To improve on present calculations, the pre- 
cision on the VP dispersion integrals is required to be 
better than 1%. More precise experimental data in the 
7r + 7r~(7) channel are needed, such that systematic un- 
certainties on the cross sections that are correlated over 
the relevant mass range are kept well below the percent 
level. 
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In this paper an analysis of the process e + e~ — > 
ir + ir~ based on data collected with the JBaBAR ex- 
periment is presented. In addition, as a cross-check of 
the analysis, we measure the e + e~ — > /i + /i~(7)7 cross 
section on the same data and compare it to the QED 
prediction. The reported results and their application to 
the 7T7T contribution to the muon magnetic anomaly have 
been already published in shorter form [7[. 



B. The ISR approach 

The initial-state radiation (ISR) method has been pro- 
posed @-[ll| as a novel way to study e+e annihilation 
processes instead of the standard point-by-point energy- 
scan measurements. The main advantage of the ISR ap- 
proach is that the final-state mass spectrum is obtained in 
a single configuration of the e + e~ storage rings and of the 
detection apparatus, thus providing a cross section mea- 
surement over a wide mass range starting at threshold. 
Consequently, a better control of the systematic errors 
can be achieved compared to the energy-scan method, 
which necessitates different experiments and colliders to 
cover the same range. The disadvantage is the reduc- 
tion of the measured cross section, which is suppressed 
by one order of a. This is offset by the availability of 
high-luminosity e + e~ storage rings, primarily designed 
as B and K factories in order to study CP violation. 

In the ISR method, the cross section for e + e~ — > X at 



the reduced energy v s' = mx , where X can be any final 
state, is deduced from a measurement of the radiative 



Xj, where the photon is emitted by the 
particle. The reduced energy is related 



cm. energy. In this analy- 
s' ranges from the two-pion 



process e + e~ - 
initial e + or e~ 
to the energy E* of the ISR photon in the e + e~ center- 
of-mass (cm.) frame by s' = s(l — 2E*/^/s) 1 where s 
is the square of the e + e 
sis, s ~ (10.58 GeV) 2 and 
production threshold to 3 GeV. Two-body ISR processes 
e + e _ — > Xj with X — 7r + 7r~(7) and X = /i + /i~(7) are 
measured, where the ISR photon is detected at large an- 
gle to the beams, and the charged particle pair can be 
accompanied by a final-state radiation (FSR) photon. 

Figure [T] shows the Feynman diagrams relevant to this 
study. The lowest-order (LO) radiated photon can be ei- 
ther from ISR or FSR. In the muon channel, ISR is dom- 
inant in the measurement range, but the LO FSR con- 
tribution needs to be subtracted using QED. In the pion 
channel, the LO FSR calculation is mo del- dependent, but 
the contribution is strongly suppressed due to the large 
s value. In both channels, interference between ISR and 
FSR amplitudes vanishes for a charge-symmetric detec- 
tor. 

In order to control the overall efficiency to high pre- 
cision, it is necessary to consider higher-order radia- 
tion. The next-to-leading-ordcr (NLO) correction in 



a amounts to about 4% [12[ with the selection used 



for this analysis, while the next-to-next-to-leading-order 
(NNLO) correction is expected to be at least one order 




FIG. 1: The generic Feynman diagrams for the processes rele- 
vant to this study with one or two real photons: lowest-order 
(LO) ISR (top left), LO FSR (top right), next-to-leading or- 
der (NLO) ISR with additional ISR (bottom left), NLO with 
additional FSR (bottom right). 



of magnitude smaller than NLO. Most of the higher- 
order contributions come from ISR and hence are inde- 
pendent of the final state. As the cross section is mea- 
sured through the mr/fj,fi ratio, as explained below, most 
higher-order radiation effects cancel and NLO is sufficient 
to reach precisions of 10~ 3 . As a result, the selection 
keeps 7T7T77 (/i/177) as well as 7T7T7 (/i/i7) final states, 
where the additional photon can be either ISR or FSR. 



C. Cross section measurement through the irn / ' 
ratio 

The cross section for the process e + e~ — > X is related 
to the \f~s' spectrum of e + e _ — > X71SR events through 



dN 



jrcff 



Ex 



(i) 



where dL^R / d\fs' is the effective ISR luminosity, £x 7 is 
the full acceptance for the event sample, and a x is the 
'bare' cross section for the process e + e~ — > X (includ- 
ing additional FSR photons), in which the leptonic and 
hadronic vacuum polarization effects are removed. 

Eq. (TTJ) applies equally to X = 7^(7) and X = ^(7) 
final states, so that the ratio of cross sections is directly 
related to the ratio of the pion to muon spectra as a func- 
tion of Specifically, the ratio i?exp(v / s 7 ) of the pro- 
duced 7T7t(7)7isr and w<(7)7isr spectra, obtained from 
the measured spectra corrected for full acceptance, can 
be expressed as: 



-Rcxp(vs') 



dN pr ° d , 

7T7T(7)7iSR 

dN pm ,\ 

W7)7lSR 



r(7) 



Ws') 



(1 



(1 



^FSrH 1 + ^add.FSR) 



(2) 

(3) 
(4) 
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The 'bare' ratio R° (no vacuum polarization, but ad- 
ditional FSR included), which enters the VP dispersion 
integrals, is given by 



CTpt(vs') 



(5) 



where cr pt = Ana 2 /3s' is the cross section for pointlike 
charged fermions. The factor (1 + ^fsr) correc ts for the 
lowest-order FSR contribution, including possibly addi- 
tional soft photons, to the e + e~ — > /i + /i~7 final state, as 
is explicitly given in Eq. (fT5)l . No such factor is included 
for pions because of the negligible LO FSR contribution 
(see Sect. HX HIT) . The factor (1 + 6%£ d FSR ) corrects for 
additional FSR in the e + e~ — > process at y/s', as 

is explicitly given in Eq. (fTT)|) . 

In this analysis, we use a procedure strictly equivalent 
to taking the ratio R exp (y/s'), namely we measure the 



• the ISR photon efficiency cancels out; 

• the vacuum polarization also cancels out. 

Furthermore the Monte Carlo generator and the detector 
simulation are only used to compute the acceptance of the 
studied X71SR processes, with X = 7^(7), /i/x(7). The 
overall systematic uncertainty on the mr cross section 
is reduced, because some individual uncertainties cancel 
between pions and muons. 



II. ANALYSIS OUTLINE 

A. The BABAR detector and data samples 

The analysis is based on 232 fb _1 of data collected with 
the BABAR detector at the SLAC PEP-II asymmetric- 



nance. The BABAR detector is described in detail else- 



c}J , s(-yV) cross section using Eq. ([T]) in which the effec- energy e + e~ storage rings operated at the T(AS) reso- 
tive ISR luminosity is obtained from the mass spectrum 
of produced /x/x(7)7isr events divided by the ^(^(Vo ) 
cross section computed with QED. The ISR luminosity 
measurement is described in detail in Sectior fVIIIFI 

This way of proceeding considerably reduces the uncer- 
tainties related to the effective ISR luminosity function 
when determined through 



Jfcff 
"-^ISR 



L, 



dW fa(s') 



(W) 



(6) 



Eq. ([6]) relies on the e + e luminosity measurement (L ee ) 
and on the theoretical radiator function dW/dy/s'. The 
latter describes the probability to radiate an ISR photon 
(with possibly additional ISR photons) so that the pro- 
duced final state (excluding ISR photons) has a mass y/s'. 
It depends on \/s, on y/s', and on the angular range (6^ in , 
6>* lax ) of the ISR photon in the e + e~ cm. system. For 
convenience, two factors that are common to the muon 
and pion channels are included in the effective luminosity 
definition of Eq. (|6]): i) the ratio of e 7lSR , the efficiency to 
detect the main ISR photon, to the same quantity eJ^^ 
in simulation, and ii) the vacuum polarization correction 
(a(s')/a(0)) 2 . The latter factor is implicitly included in 
the effective luminosity deduced from ////(7)tisr data us- 
ing Eq. (JT|), while the former, which cancels out in the irir 
to fifi ratio, is ignored in Eq. ([T]). As an important cross- 
check of the analysis, hereafter called the QED test, we 
use Eq. (|6]), together with Eq. ([T]), to measure the muon 
cross section and compare it to the QED prediction. 

Many advantages follow from taking the R e ^ P (y/s') ra- 
tio: 

• the result is independent of the BABAR luminosity 
L ee measurement; 

• the determination of the ISR luminosity comes 
from the muon data, independently of the num- 
ber of additional ISR photons, and thus does not 
depend on a theoretical calculation; 



where 13j. Charged-particle tracks are measured with a 
five-layer double-sided silicon vertex tracker (SVT) to- 
gether with a 40-layer drift chamber (DCH) inside a 
1 .5 T superconducting solenoid magnet. Photons are as- 
sumed to originate from the primary vertex defined by 
the charged tracks of the event and their energy is mea- 
sured in a CsI(Tl) electromagnetic calorimeter (EMC). 
Charged-particle identification (PID) uses the ionization 
losses dE/dx in the SVT and DCH, the Cherenkov ra- 
diation detected in a ring-imaging device (DIRC), the 
shower energy deposit in the EMC (-E C al) and the shower 
shape in the instrumented flux return (IFR) of the mag- 
net. The IFR system is constructed from modules of 
resistive plate chambers interspaced with iron slabs, ar- 
ranged in a configuration with a barrel and two endcaps. 



B. Monte Carlo generators and simulation 

Signal and background ISR processes e + e~ — > X71SR 
are simulated with a Monte Carlo (MC) event gener- 
ator called AfkQed, which is based on the formalism 
of Ref. [f4|. The main ISR photon, 7isr, is generated 
within the angular range [9* min = 20°, 6> t * nax = 160°] in 
the cm. system, bracketing the photon detection range 
with a margin to account for finite resolution. Addi- 
tional ISR photons are generated with the structure func- 
tion method [f5[ , and additional FSR photons with PHO- 
TOS [n| . Additional ISR photons are emitted along the 
e + or e~ beam particle direction. A minimum mass 
m x 7isr > 8 GeV/ c 2 is imposed at generation, which 
places an upper bound on the additional ISR photon en- 
ergy. Samples corresponding to 5 to 10 times the num- 
ber of data events are generated for the signal channels. 
The more accurate Phokhara generator [17| is used at 
the 4-vector level to study some effects related to addi- 
tional ISR photons. Background processes e + e~ — > qq 
(q = u,d,s,c) are generated with JETSET [llj], and 
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e+e~ t+t~ with KORALB [lj|. The response of the 
BAB4R detector is simulated with GEANT4 [20| . 



C. Analysis method 

The 7T7t(7)7isr and /z/u(7)7isr processes are measured 
independently with full internal checks and the ratio 
-Rcx P (v / s 7 ), which yields the measured a s ^ v ^{-/s i ) cross 
section, is only examined after these checks are success- 
fully passed. One of the most demanding tests is the ab- 
solute comparison of the /UA*(7)7isr cross section, which 
uses the BABAR L ee luminosity, with the NLO QED pre- 
diction (QED test). 

After preliminary results were presented from the blind 
analysis [21[ , a few aspects of the analysis were revisited 
to refine some effects that had been initially overlooked, 
mostly affecting the correlated loss of muon identifica- 
tion for both tracks. While the final measurement is not 
a strictly blind analysis, all studies are again made inde- 
pendently for muons and pions and combined at the very 
end. 

The selected events correspond to a final state with 
two tracks and the ISR candidate, all within the detec- 
tor acceptance, as described in Section Hill Kinematic fits 
provide discrimination of the channels under study from 
other processes. However the separation between the dif- 
ferent two-prong final states (including K + K~ (7)tisr) 
relies exclusively on the identification of the charged par- 
ticles. Thus particle identification plays a major role 
in the analysis. This is the subject of Section IIVDI 
Background reduction and control of the remaining back- 
ground contributions are another challenge of the anal- 
ysis, in particular in the pion channel away from the p 
resonance. This is discussed in Section fVTl 

The determination of the y/s' spectrum is described in 
Section IVIIl The relevant final-state mass is (m^) 
when there is additional ISR or no additional radiation, 
or m 7r7r7 (m MM7 ) in the case of additional FSR. The y/s* 
spectrum is obtained from the observed (m w ) dis- 
tributions through unfolding (Sect. EH A[) . 

Although selection of the final state of two-body ISR 
processes is rather simple, the main difficulty of the anal- 
ysis resides in the full control of all involved efficiencies. 
Relying on the simulation alone cannot provide the re- 
quired precision. The simulation is used in a first step 
in order to incorporate in a consistent way all effects en- 
tering the final event acceptance. Corrections for data- 
to-MC differences are obtained for each efficiency using 
dedicated studies performed on the data and simulation 
samples. The main contributions for these corrections 
originate from trigger, tracking, particle identification, 
and the x 2 selection of the kinematic fits, so that the 
corrected efficiency is 

( £ data\ / data \ / £ data\ / £ data \ 
tIrT TW 1 ) Jf ) 3Tc • ( ? ) 
fc trig / \ fc track/ V fc PID / \ fc X 2 / 



The corrections Ci = ( J are reviewed in turn in the 

following sections (Sect . HVl and |V|) . They are applied as 
mass-dependent corrections to the MC efficiency. They 
amount to at most a few percent and are known to a 
few permii level or better. Efficiency measurements are 
designed to avoid correlations between the Cj. Further 
data-to-MC corrections deal with second-order effects re- 
lated to the description of additional ISR in the genera- 
tor, which was found inadequate at the level of precision 
required for this analysis. 



III. EVENT SELECTION 



A. Topological selection 

Two-body ISR events are selected by requiring a pho- 
ton candidate with E* > 3 GeV and laboratory polar 
angle in the range 0.35 — 2.4 rad, and exactly two tracks 
of opposite charge, each with momentum p > 1 GeV/c 1 
and within the angular range 0.40 — 2.45 rad. A photon 
candidate is defined as a cluster in the EMC, with energy 
larger than 0.02 GeV, not associated to a charged track. 
If several photons are detected, the main ISR photon is 
assumed to be that with the highest E*; this results in 
an incorrectly assigned ISR photon in less than ~ 10~ 4 
of the events, mostly due to the ISR photon loss in inac- 
tive areas of the EMC. The track momentum requirement 
is dictated by the fall-off of the muon-identification effi- 
ciency at low momenta. The tracks are required to have 
at least 15 hits in the DCH, and originate within 5 mm of 
the collision axis (distance of closest approach doca xy < 
5 mm) and within 6 cm from the beam spot along the 
beam direction (|A Z | < 6 cm). They are required to ex- 
trapolate to the DIRC active area, whose length further 
restricts the minimum track polar angle to ~ 0.45 rad. 
Tracks are also required to extrapolate to the IFR ac- 
tive areas that exclude low-efficiency regions. An addi- 
tional veto based on a combination of E ca \ and dB/dx, 
((S cal /p-l)/0.15) 2 + ((dB/da;DCH-690)/150) 2 < 1), 
reduces electron contamination. Events can be accom- 
panied by any number of 'bad' tracks, not satisfying the 
above criteria, and any number of additional photons. 
To ensure a rough momentum balance at the preselection 
level (hereafter called 'preselection cut'), the ISR photon 
is required to lie within 0.3 rad of the missing momentum 
of the tracks (or of tracks plus other photons). 



Unless otherwise stated, starred quantities are measured in the 
e + e~ cm. and un-starred quantities in the laboratory. 
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FIG. 2: (color online). The 2D-* 2 distribution for 
7T7r(7)7isR (data) for 0.5 < TO W7r < 1.0GeV/c 2 , where dif- 
ferent interesting regions are defined. 



FIG. 3: (color online). The 2D-x 2 distribution for 
A t A*(7)7iSR (data) for 0.5 < m MM < 1.0GeV/c 2 , where the 
signal and background regions are indicated. 



B. Kinematic Fit description and \ 2 selection 

For both the /i/17 and irn^/ processes, the event defini- 
tion is enlarged to include the radiation of one photon in 
addition to the already-required ISR photon. Two types 
of fits are considered, according to the following situa- 
tions: 

• The additional photon is detected in the EMC, in 
which case its energy and angles can be readily 
used in the fit: we call this a 3-constraint (3C) 
FSR fit, although the extra photon can be either 
from FSR or from ISR at large angle to the beams. 
The threshold for the additional photon is kept low 
(20MeV). This can introduce some background, 
but with little effect as the fit in that case would 
not be different in practice from a standard fit to 
the a*m(7)7isr (tttt^Tisr) hypothesis. 

• The additional photon is assumed to be from ISR 
at a small angle to the beams. Since further infor- 
mation 2 is not available, it is presumed that the 
extra photon is perfectly aligned with either the 
e + or the e~ beam. The corresponding so-called 
2C ISR fit ignores additional photons measured in 



2 This is not strictly true as the missing photon could be com- 
pletely reconstructed if the ISR photon energy were used in the 
kinematic fit. However tests have shown that the relative quality 
of this new information does not permit a significant improve- 
ment for the fitted direction of the additional ISR photon over 
the collinear assumption. 



the EMC and determines the energy of the fitted 
collinear ISR photon. 

In both cases the constrained fit procedure uses the 
ISR photon direction and the measured momenta and 
angles of the two tracks with their covariance matrix in 
order to solve the four energy-momentum conservation 
equations. The measured energy of the primary ISR pho- 
ton is not used in either fit, as it adds little information 
for the relatively low masses involved. 

Each event is characterized by two x 2 values, Xfsr 
and Xisr from the FSR and ISR fits, respectively, which 
are examined on a 2D-plot. In practice the quantities 
ln(x 2 + 1) are used so that the long tails can be prop- 
erly visualized (Figs. [3J Events without any extra 
measured photons have only the Xisr value and they are 
plotted separately on a line above the Xfsr overflow. In 
case several extra photons are detected, FSR fits are per- 
formed using each photon in turn and the fit with the 
smallest Xfsr 1S retained. The muon (pion) mass is as- 
sumed for the two charged particles, according to the 
selected channel, and in the following studies and final 
distributions, the /Lt/i (7T7t) mass is obtained using the fit- 
ted parameters of the two charged particles from the ISR 
fit if Xisr < Xfsr ano - fr° m the FSR fit in the reverse 
case. 

It is easy to visualize the different interesting regions 
in the 2D-x 2 plane, as illustrated in Fig. [5] for 7T7r(7)7isR 
data. Most of the events peak at small values of both x 2 , 
but the tails along the axes clearly indicate events with 
additional radiation: small-angle ISR along the Xfsr ax ^ s 
(with large ISR energies at large values of Xfsr)' or FSR 
or large-angle ISR along the Xisr ax i s (with large addi- 
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tional radiation energies at large values of Xisr ) ■ Events 
along the diagonal do not satisfy either hypothesis and 
result from resolution effects for the pion tracks (also 
secondary interactions) or the primary ISR photon, or 
possibly additional radiation of more than one photon. 
Multibody background populates the region where both 
X 2 are large and consequently a background region is de- 
fined in the 2D-% 2 plane. This region is optimized as a 
compromise between efficiency and background contami- 
nation in the signal sample, aiming at best control of the 
corresponding systematic uncertainties. 

The x 2 criteria used in the pion analysis depend on the 
7T7T mass region considered. The m-n-n region between 0.5 
and 1 GeV/c 2 is dominated by the p resonance. The cor- 
responding large cross section provides a dominant con- 
tribution to vacuum-polarization dispersion integrals, so 
it has to be known with small systematic uncertainties. 
Also background is expected to be at a small level in 
this region. These two considerations argue for large ef- 
ficiencies, in order to keep systematic uncertainties suffi- 
ciently low. Therefore a loose x 2 criterion is used, where 
the physical (accepted) region corresponds to the left of 
the contour outlined in Fig. [2j excluding the BG-labeled 
region. The same loose x 2 criterion is applied for the 
M/«(7)7iSR analysis (Fig.©. 

The pion form factor decreases rapidly away from the 
p peak, while the backgrounds vary slowly with the ttt: 
mass. The multihadronic background in the physical 
sample becomes excessively large if the x 2 criterion as 
used in the p region is applied, and it is necessary to 
tighten the selection of 7T7t(7)7isr events. Figure 2] shows 
the tight x 2 selection boundary ln(x 2 SR + l) < 3 chosen to 
reduce multihadronic background, and the 2D-% 2 distri- 
butions for masses below and above the central p region. 
The tight x 2 criterion retains events with additional ISR 
since this region in the % 2 plane is free of multihadronic 
background. The reduced efficiency on signal from the 
tight selection results in a larger relative uncertainty, but 
this is still acceptable considering the much smaller con- 
tribution from the p tails to the dispersion integral. 

Efficiencies and systematic uncertainties resulting from 
the loose and tight x 2 selection criteria are discussed in 
Section EE 



IV. EFFICIENCY STUDIES (I) 

To achieve the required precision for the cross section 
measurement, efficiencies are validated with data at ev- 
ery step of the event processing, and mass-dependent 
data/MC corrections are determined. This necessitates 
specific studies on data control samples whose selection 
criteria are designed to minimize biases on efficiency mea- 
surements. Residual effects are estimated and included 
in the systematic errors. 



A. Efficiency-dedicated event selection and 
kinematic fit 

For trigger and tracking efficiency studies, a dedicated 
selection of /j + /i~7isr and 7t + 7t~7isr events is devised 
that only requires one reconstructed track (called 'pri- 
mary'), identified as a muon or pion, and the ISR pho- 
ton. A 1C kinematic fit is performed and the momen- 
tum vector of the second muon(pion) is predicted from 
4-momentum conservation. Standard track selection is 
applied to the primary track and the predicted track is 
required to be in the acceptance. 



B. Trigger and filtering 

A number of trigger conditions are imposed at the 
hardware (LI) and online software (L3) levels, as well as 
in a final filtering, before an event is fully reconstructed 
and stored in the BABAR data sample. They are common 
to all BABAR analyses, and hence are not specifically de- 
signed to select ISR events. Since individual trigger and 
filter line responses are stored for every recorded event, 
efficiencies can be computed by comparing the response 
of trigger lines, after choosing lines that are as orthog- 
onal and as efficient as possible. Trigger efficiencies are 
determined on data and simulation samples, after ap- 
plying identical event selections and measurement meth- 
ods, and data/MC corrections C tr i g are computed from 
the comparison of measured efficiencies on background- 
subtracted data and signal MC. Once the physics ori- 
gins of inefficiencies are identified, uncertainties are es- 
timated through studies of biases and data-to-MC com- 
parison of distributions of relevant quantities. Efficien- 
cies and data/MC corrections are measured separately 
for the pion and muon channels. 

Trigger efficiencies are determined on samples unbi- 
ased with respect to the number of tracks actually recon- 
structed, to avoid correlations between trigger and track- 
ing efficiency measurements. In practice, one- and two- 
track samples are sufficient and consequently the trig- 
ger control samples are selected through the dedicated 
1C kinematic fit described above. Because of the loose 
requirement with respect to tracking, the data samples 
contain backgrounds with potentially different trigger ef- 
ficiencies to that of the signal. These backgrounds are 
studied with simulation and are then subtracted. To 
obtain data samples that are as pure as possible, crite- 
ria tighter than the standard track selection are applied 
to the primary track, including tight PID identification. 
Possible biases resulting from the tighter selection are 
studied and accounted for in the systematic errors. Back- 
ground contributions are subtracted from the data spec- 
tra using properly-normalized simulated samples, and, if 
necessary, with data/MC correction of the trigger effi- 
ciencies in an iterative procedure. 

The data/MC corrections for the LI trigger are found 
to be at a few x 10~ 4 level for muon and pion events. 
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FIG. 4: (color online). The 2D-x 2 distributions in 7T7t(7)7isr data: (left) below the central p region (m nn < 0.5 GeV/c 2 ); (right) 
above the central p region (1. < m vv < 2. GeV/c 2 ). The line indicates the boundary for the tight x 2 selection. 



The L3 level involves a track trigger (at least one track 
is required) and a calorimetric trigger (demanding at 
least one high-energy cluster and one low-energy clus- 
ter). Both of them are efficient for 7T7T7isr events. For 
A*/^7isr events, the small efficiency of the calorimetric 
trigger limits the statistical precision of the track-trigger 
and overall efficiency measurements. Furthermore, a cor- 
related change of the two trigger line responses for close- 
by tracks induces both a non-uniformity in the efficiency 
and a bias in the efficiency measurement. This origi- 
nates from the overlap of tracks in the drift chamber and 
of showers in the EMC, which induces a simultaneous 
decrease in the track-trigger efficiency and an increase 
in the calorimetric-trigger efficiency. Overlap is a major 
source of overall inefficiency and difference between data 
and simulation, necessitating specific studies. The cor- 
rection to the MC L3 trigger efficiency is small for pions, 
about 2 x 10 -3 at the p peak, and known to a precision 
better than 10~ 3 . The data/MC correction Ct r i g is larger 
in the /i/x(7)7isR channel, due to the dominant role of the 
track trigger, about 1% at a pip, mass of 0.7GeV/c 2 , and 
known to a precision of 3 x 10 -3 (Fig. [S]top). Uncertain- 
ties, which increase to 5 x 10~ 3 at the maximum overlap 
("V^i ~ 0.4GeV/c 2 ), are mostly statistical in nature. 

The offline event filtering involves a large number 
of specific selections, including dedicated ppj filters. 
Whereas the inefficiency and its correction are negligi- 
ble for muons, some inefficiency at the filtering stage 
is observed for pion events, mostly at low m^ir mass. 
This originates again from the overlap of tracks in the 
DCH and hadronic showers in the EMC. The correction 
Ctrig to the 7T7t(7)7isr cross section (Fig. [SJ bottom) is 
found to be (1.0 ± 0.3 ± 0.3)% at ~ 0.4GeV/c 2 and 



(2.9 ± 0.1 ± 1.0) x 10" 3 at the p peak, where the first 
error is statistical and the second systematic. Beyond 
1.5GeV/c 2 , the background level precludes a significant 
measurement of the efficiency with data and no correc- 
tion is applied; a systematic error of 0.4 x 10~ 3 is assigned 
in the high mass range, equal to the inefficiency observed 
in MC. Due to imperfect simulation of hadronic showers, 
filtering is the major source of trigger systematic uncer- 
tainties in the pion channel. 

Systematic errors due to the trigger and filter are re- 
ported in Tables |TT] and |Vl for muon and pion channels, 
respectively. 



C. Tracking 

The tracking control samples of ptpijisR (""""7isr) 
events are selected through the efficiency-dedicated 1C 
fit described above. The rate of predicted tracks that 
are actually reconstructed in the tracking system, with a 
charge opposite to that of the primary track, yields the 
tracking efficiency. 

To ensure the validity of the measurement, further cri- 
teria are applied to the tracking sample in addition to 
the kinematic fit. To enhance purity, a 7r° veto is applied 
if a pair of additional photons in the event can form a 
7T° candidate with mass within 15 MeV/c 2 of the nominal 
mass. This tt° veto is not applied to the pion tracking 
sample, because of the bias it would introduce on the in- 
efficiency related to secondary interactions. The events 
are required to pass the triggers and online filter and 
are selected without specific requirements on the second 
reconstructed track, if any. Biases affecting the tracking- 
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FIG. 5: The data/MC event trigger and filter correction Ctrig 
for the pp(y)yisR (top) and 7T7t(7)7isr (bottom) cross sections 
as a function of the fifi and rnr masses, respectively. Statistical 
errors only. 



efficiency measurement introduced by the primary-track 
selection or event-level background-rejection criteria are 
identified and studied with simulation, and evaluated 
with data. 

The predicted track is required to lie within the track- 
ing acceptance, taking into account the effect of angular 
and momentum resolution. The method therefore deter- 
mines the efficiency to reconstruct a given track in the 
SVT+DCH system within a specified geometrical accep- 
tance, no matter how close or distant this track is with 
respect to the expected one, due for instance to decays or 
secondary interactions. However, the possible mismatch 
in momentum and/or angles affects the full kinematic 



reconstruction of the event, and its effect is included in 
the efficiency of the corresponding % 2 selection applied to 
the physics sample (Sect. IVB4p . Likewise, effects from 
pion decays in the detector volume are included in the 
pion-identification efficiency. 

The individual track efficiency is determined assuming 
that the efficiencies of the two tracks are uncorrelated. 
However, the tracking efficiency is observed to be sharply 
reduced for overlapping tracks in the DCH, as measured 
by the two-track opening angle (A(f>) in the plane trans- 
verse to the beams. Not only is the individual track ef- 
ficiency locally reduced, but a correlated loss of the two 
tracks is observed. In addition, as the final physics sam- 
ple is required to have two and only two good tracks 
with opposite charge, the understanding of the tracking 
involves not only track losses, but also the probability to 
reconstruct extra tracks as a result of secondary interac- 
tions with the detector material or the presence of beam- 
background tracks. The full tracking efficiency is then 
the product of the square of the single-track efficiencies, 
the probability for not losing the two tracks in a corre- 
lated way (loss probability = /o), and the probability for 
not having an extra reconstructed track (loss probability 
= /a). The event correction Ct ra ck to be applied to the 
MC is the corresponding product of the data/MC ratios 
of each term. The mass-dependent quantities fo and 
are in the (0.5 — 2.5) x 10 -3 range. 

For muons the single-track inefficiency and the 
data/MC correction are driven by the DCH overlap ef- 
fect. At the maximum overlap (rn w ~ 0.4GeV/c 2 ) 
the inefficiency reaches 1.7% in simulation, but 2.5% in 
data, while the intrinsic reconstruction, measured for 
non-overlapping tracks, accounts for an inefficiency of 
2.5 x 10~ 3 in data, and 5 x 10~ 4 in simulation. 

Because of backgrounds, the pion tracking efficiency 
can be obtained directly from data only in the p peak re- 
gion, from 0.6 to 0.9GeV/c 2 . The main sources of track 
loss are identified: the track overlap in the DCH and the 
secondary interactions. The two effects are separated us- 
ing the A<fi distribution. This two-component model is 
used to extrapolate the inefficiency to mass regions out- 
side the p peak. Results for pions are qualitatively similar 
to those for muons, with inefficiencies driven by the track 
overlap effects. The intrinsic track inefficiency is domi- 
nated by secondary interactions (2.2% in data and 1.7% 
in simulation) and is thus much larger than for muons. 
Near 0.4 GeV/c 2 the track inefficiency is determined to be 
6.2% in data and 4.7% in simulation. Above 1.2GeV/c 2 
for pions, where the data/MC correction is not expected 
to vary significantly, a systematic uncertainty of about 
0.3% is assigned. 

The final corrections Ct ra ck to the /i/i(7)7isR and 
7T7t(7)7isr cross sections are presented in Fig. O Ct rac k 
differs from unity by about 1.6% (3.0%) at 0.4GeV/c 2 , 
and by 0.8% (1.5%) at 1 GeV/c 2 for muons (pions). Sta- 
tistical uncertainties from the efficiency measurements 
are indicated by point-to-point errors. Systematic un- 
certainties are estimated from the study of biases in the 
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method, determined using the simulation and calibrated 
with data-to-MC comparison of distributions character- 
istic of the physics source of the bias. Systematic un- 
certainties amount to 0.8 x 10~ 3 for muons in the mass 
range from 0.4 to 1.0GeV/c 2 , and are about a factor of 
two larger outside. For pions, the systematic uncertainty 
of the correction is 1.1 x 10 -3 in the 0.6 — 0.9 GeV/c 2 mass 
range, increasing to 2.1 x 10~ 3 (0.4 - 0.6GeV/c 2 ), 3.8 x 
10" 3 (below 0.4GeV/c 2 ), 1.7 x 10~ 3 (0.9- 1.2GeV/c 2 ), 
and 3.1 x 10~ 3 (above 1.2GeV/c 2 ). 
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D. Particle identification 

The method to determine the PID efficiencies makes 
use of the x + x~yisu sample itself, where one of the pro- 
duced charged x particles (x = /i, ir, K) is tagged using 
strict identification criteria, and the second ('opposite') 
track identification is probed ('tag-and-probe' method). 
The events are selected through a 1C kinematic fit that 
uses only the two tracks, with assumed mass m x . The 
requirement Xxx < 15 is applied to reduce multihadronic 
background. In this way the ensemble of opposite tracks 
constitutes a pure x sample to be subjected to the iden- 
tification process. The residual small impurity in the 
data samples is measured and corrected in the efficiency 
determination. The same analysis is performed on MC 
samples of pure x + x~^isk events, and data/MC correc- 
tions Cpid are determined for each x type, as explained 
below. Since the PID efficiency measurement relies on 
two-track events that have passed the triggers, Cpid is 
not correlated with Ct r i g or Ct rac k, as required by Eq. J?}. 



1. Particle identification classes 



Particle ID measurements in this analysis aim to ob- 
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FIG. 6: The data/MC event tracking correction Ctrack for the 
wMtisr (top) and 7T7t(7)7isr (bottom) cross sections as a 
function of the nfi and tttv masses, respectively. 



tain from data the values for all the elements e. 



of 



the efficiency matrix, where i is the true e, /i, 7r, or K 
identity and 'j' is the assigned ID from the PID proce- 
dure (Table Protons (antiprotons) are not included in 
the particle hypotheses because the pp^f final state oc- 
curs only at a very small rate [22|. This contribution 
is estimated from simulation, normalized to data, and 
subtracted statistically from the mass spectra. 

We identify muon candidates by applying criteria on 
several discriminant variables related to the track, such 
as the energy deposition E ca \ in the EMC, and the track 
length, hit multiplicity, matching between hits and ex- 
trapolated track in the IFR. This defines the hyq selector. 
The if id selector is constructed from a likelihood func- 
tion using the distributions of cYE/dx in the DCH and of 
the Cherenkov angle in the DIRC. The electron identifi- 
cation relies on a simple E ca \/p > 0.8 requirement. As 
most of the electrons are vetoed at the preselection level, 
their fraction in the pion sample is generally small. Their 
contribution is completely negligible in the muon sample. 

In addition to physical particle types, we assign an 
ID type of '0' if the number of DIRC photons associ- 
ated with the track (./Vdirc) is insufficient to define a 
Cherenkov ring, thus preventing tt-K separation. The 
ID classes defined in Table |T] constitute a complete and 
orthogonal set that is convenient for studying cross-feed 
between different two-prong ISR final states. 

The 7T-ID selector is a set of negative conditions since 
no set of positive pion-ID criteria was found that provided 
both sufficient efficiency and purity. Pion candidates are 
tracks that do not satisfy any of the other ID class re- 
quirements. In this sense the pion-ID is sensitive to the 



14 



TABLE I: Definition of particle ID types (first column) us- 
ing combinations of experimental conditions (first row) : "+" 
means "condition satisfied" , "— " means "condition not satis- 
fied" , an empty box means "condition not applied" . The con- 
ditions jUm and Km correspond to cut-based and likelihood- 
based selectors, respectively. The variables -/Vdirc and E ca \ 
correspond to the number of photons in the DIRC and the 
energy deposit in the EMC associated to the track, respec- 
tively. 

1 | mid | ffcai/p > 0-8 1 jVpiRc < 2 1 Km | 
I y I + I | I I 

| '0' | - I - I + I I 

l2g] - I - I - I + I 



problems affecting the identification of all the other par- 
ticle types. 



2. 'Hard' pion-ID definition 

The standard 7T-ID definition in Table U is part of a 
complete set of exclusive PID conditions that is conve- 
nient when backgrounds in the pion candidate sample are 
expected to be manageable. However, in some cases the 
standard 7T-ID algorithm does not deliver sufficient pion 
purity. One such case concerns the purity of the tagged 
pion in the tag-and-probe pion pair, which is crucial 
in the determination of 7T-ID efficiencies (Sect. IIVD 5p . 
Improved pion purity is also necessary to reduce back- 
grounds for the 7T7r cross section measurement in mass 
regions in the tails of the p peak. 

A tighter pion-ID selector is thus developed, which we 
call the 'hard pion' (iTh) selector, to improve the rejection 
of muons and electrons that are misidentified as pions 
with the standard definition. The n^-ID is based on two 
likelihood functions P^/^ and P T / e : P-k/^j, uses the EMC 
deposit E ca \ associated with the track and the penetra- 
tion of the track into the IFR, while P^/ e uses E ca \ and 
the measurements of (dE/dx)DCH and (dE/da;)svT as 
a function of momentum. Tracks with likelihoods close 
to correspond to pions while P-n/n ~ 1 (P-x/e ~ 1) are 
muon-like (electron- like) . Reference distributions used in 
the likelihoods are obtained from simulation, with cor- 
rections determined from data control samples in mass 
ranges that ensure backgrounds are negligible. The pure 
muon sample used for PID efficiency (see below), and 
the sample identified as 'e7r' in the very high mass range 
(m„„ > 5 GeV/c 2 ), provide the reference distributions for 
muons and electrons, respectively. The 7777(7)71311 sam- 
ple, with m^T,- restricted to the p peak and with both pi- 
ons satisfying the standard 7T-ID, provides the reference 
distributions for pions. 



3. PID measurements with the muon samples 

The method used to determine the muon ID efficiency 
utilizes of the p,[i (7)7isr sample itself, where one of the 
tracks is tagged as a muon using the /iin selector defined 
above, and the opposite track is probed. The sample is 
restricted to m MM > 2.5 GeV/c 2 to reduce the non-/^ back- 
ground to the (l.lztO.l) x 10 -3 level, so that the ensemble 
of opposite tracks constitutes a pure muon sample. 

The IFR performance at the time the data for this 
paper were collected was non-uniform across the detec- 
tor and deteriorated with time 3 . In order to map the 
PID efficiency, the track to be probed is extrapolated to 
the IFR. Local coordinates (vi,i>2) of the impact point 
are defined depending on the IFR geometry (barrel or 
endcaps). Efficiency maps are obtained for each of the 
four data-taking periods used in the analysis. The gran- 
ularity of the 3-dimensional (3D) maps is optimized as a 
function of momentum and local coordinates (p, v±,V2), 
so that local variations of efficiencies are described with 
significant statistical precision. 

The low-efficiency regions in the IFR are removed in 
order to keep as active areas only the regions where the 
Pid efficiency was reasonably homogeneous. Removed 
portions include the crack areas between modules and 
some parts of the nominal active region where the IFR 
performance was strongly degraded. The definition of 
removed regions is run-dependent: in the first running 
period only cracks are removed (about 13% of the IFR 
solid angle), while in the fourth period an additional ~ 
15% is eliminated. 

Due to the mass restriction applied to the muon con- 
trol sample, the 3D-maps provide the identification effi- 
ciency for isolated muon tracks. They parametrize the 
local performance of the IFR at the track impact point. 
However, at \i\x masses less than 2.5GeV/c 2 , tracks can 
become geometrically close to each other within the IFR 
and their respective ID efficiencies can be significantly af- 
fected. First, the efficiency is reduced with respect to the 
isolated track efficiency because the combination of the 
two sets of hits causes some of the criteria that enter the 
Pyq selector to fail. Also the single- hit readout of the two- 
dimensional strip structure of the IFR chambers leads to 
losses. Second, track overlap leads to a correlated loss of 
PID for both tracks, not accounted for by the product of 
their uncorrelated single-track inefficiencies registered in 
the maps. 

The loss of efficiency and the correlated loss effects 
are studied and evaluated in data and in simulation us- 
ing the two-track physics sample. Since the pion back- 
ground in the data sample is large in the p mass range, 
the efficiencies are measured directly only in the mass 
regions in the resonance tails, and are then extrapolated 



3 This problem was remedied for data collected subsequently, 
through IFR detector upgrades. 
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to the p peak region (0.6 — 0.9GeV/c 2 ). Possible bias 
from this procedure is studied with simulation and a sys- 
tematic uncertainty of 2.2 x 10~ 3 is assigned. The ef- 
ficiency loss (compared to the isolated muon efficiency) 
is determined using a muon-ID tagged track as for the 
high-mass sample. Results are stored in mass-dependent 
2D-maps as a function of the differences (Av\, Ai^) be- 
tween the impact points of the two tracks in the local IFR 
coordinate system. Background from pions and kaons is 
subtracted using simulation with data/MC corrections 
for the mis-ID probabilities. The efficiency loss is max- 
imal for m M(1 ~ 0.7GeV/c 2 with a reduction of 8.4% of 
the single-track efficiency in data and 4.8% in simula- 
tion. The resulting muon-ID inefficiencies (1 — e^-^) 
measured in data and simulation at 0.4 (1.0) GeV/c 2 are 
7.7% (6.6%) and 4.2% (3.5%), respectively. 

The two-track sample with no identified muon is used 
to measure the correlated efficiency loss. Since the pion 
background is overwhelming large in that sample, the 
small di-muon component is extracted by applying the 
likelihood estimator described in Section TlV D 21 to both 
tracks. The correlated muon-ID loss is found to be 1.4% 
in data and 0.3% in simulation, at m w ~ 0.7GeV/c 2 . A 
systematic uncertainty of 1.5 x 10~ 3 is assigned to the 
data/MC correction. Both efficiency loss and correlated 
loss decrease for higher masses and vanish at 2.5 GeV/c 2 . 

The event data/MC corrections resulting from requir- 
ing muon-ID for both tracks are obtained separately for 
the different running periods. The overall correction is 
given in Fig. [7] The plotted errors are statistical only. 

Systematic errors are estimated for the different data- 
taking periods. The overall systematic error from muon- 
ID on the ai/x(7)7isr cross section amounts to 3.3 x 10~ 3 . 
Muon-ID is the largest source of uncertainty in the muon 
analysis. The dominant contribution arises from the pro- 
cedures used to estimate the efficiency loss and the cor- 
related loss of the two muons. 

The pure muon sample of opposite tracks is used to 
measure the mis-ID probabilities. The largest one is 
£ ^->'7r', very close to 1 — e^-^'. Therefore, the preceding 
results for muon-ID efficiency and data/MC corrections, 
including those for close tracks, translate to the muon 
mis-ID to 'it'. The small mis-ID probabilities into par- 
ticle types other than V are estimated from additional 
studies. The €^>k' mis-ID is smaller than 0.1% for mo- 
menta below 3 GeV/c, with a steep increase for larger mo- 
menta, reaching about 1% at 5 GeV/c, with a data/MC 
correction of 12%. Mis-ID probabilities to e^-j-'o' and 
e^i-s-'e' are 0.4% and less than 0.1%, respectively. 



4- PID measurements with the kaon samples 

For the kaon efficiency and mis-ID measurement, the 
same tag-and-probc method is used as described for the 
muons, this time with a primary track satisfying the Kid 
condition. In addition to the restriction Xkk < ^ a P~ 
plied to reduce multihadronic background, a requirement 
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FIG. 7: The data/MC event correction Cpid for muon-ID 
efficiency as a function of m M/J . 

Xkk < xi-w i s applied to reduce the pion contamination. 
The purity of the kaon sample is further enhanced by a 
restriction on the fitted ttikk mass, which must be in the 
4> window 1.01 — 1.04GeV/c 2 . The electron background 
from photon conversions in the process e + e~ — > 77, 
which populate the ttikk threshold region, is eliminated 
by a requirement on the distance (in the transverse plane) 
between the vertex of the two tracks and the beam axis. 
The purity achieved is (99.0 ±0.1)%, determined from a 
fit of the rriKK distribution in data, with cf> signal and 
background shapes taken from MC. 

The data/MC corrections for the Kyo efficiency are 
obtained as a function of track momentum. The restric- 
tion to the 4> sample imposes kinematic restrictions on the 
kaon momentum, with a lack of statistics below 1.5 GeV/c 
and above 5 GeV/c. This necessitates an extrapolation, 
which is achieved through a fit of the kinematically avail- 
able data. A sampling of the momentum-dependent cor- 
rection is performed using the K i^7isR MC simulation, 
in order to determine the event correction as a function 
of the rriKK invariant mass. 

The mis-ID probabilities depend on momentum, espe- 
cially for the K — ¥ V mis-ID, which increases strongly 
for large momenta, where the K\d selector becomes ineffi- 
cient. At 4 GeV/c the values in data are 5.8% for K — » '//', 
16.1% for K ->■ V, and 0.7% for K -> 'e'. The corre- 
sponding data/MC corrections are 0.61±0.05, 0.87±0.04, 
and 2.7 ±0.8. 



5. PID measurements with the pion samples 

The tag-and-probe method is again applied to con- 
struct a pure pion sample used to measure the pion ID 
efficiency and misidentification probabilities. To reduce 
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the backgrounds, the mass range of selected 7T7t(7)7isr 
candidates is restricted to the p peak, 0.6 < < 
0.9GeV/c 2 . To ensure the validity of the pion ID effi- 
ciency measurement, the purity of the pion ID sample is 
further increased by requiring the primary track to sat- 
isfy the 'hard pion' tag. In the restricted mass range, 
the sum of p, K and e backgrounds is reduced to the 
(3.7 ± 0.5) x 10~ 3 level. 

While for muons it is possible to measure the ID effi- 
ciencies for isolated tracks using events with a large 
mass, for pions we use events in the p region. Tracks 
in this region may overlap in one detector or another. 
Thus the 7T-ID efficiencies and mis-ID probabilities con- 
tain some average of overlap effects, which are not possi- 
ble to sort out in detail. However these effects arc much 
reduced for pions compared to muons, since showering in 
the IFR is sufficient to distinguish hadrons from muons 
and the overlap of pion showers does not degrade the 
pion-ID efficiency. 

All 7T-ID efficiencies and mis- ID probabilities are stored 
in 2D-maps as a function of the momentum and local z 
coordinate of the track extrapolated to the most rele- 
vant detector (IFR or DIRC). Biases from primary pion 
tagging and correlated two-track pion-ID loss are stud- 
ied with simulation, and verified with data in the most 
critical cases. Both effects are at the 10~ 3 level. The 
7T-ID maps are sampled to build the full event efficiency 
distributions as a function of the m nl7 mass, in data and 
MC. The event 7T-ID efficiency is weakly mass-dependent 
with typical values in data of 77.8%, 75.3%, and 77.0%, 
at masses of 0.35GeV/c 2 , 0.6GeV/c 2 , and lGeV/c 2 , re- 
spectively. 

The data/MC correction to the full event 7T-ID effi- 
ciency is shown in Fig. [8] The correction is smaller than 
the corresponding factor for muons, which reflects a lesser 
sensitivity of the 7T-ID efficiency to the IFR conditions. 
Although it has been obtained using maps determined in 
the p region, it shows only a few percent variation with 
mass, consistent with the fact that correlated ID losses 
are small. 

The systematic errors on the efficiencies come from 
the limited granularity of the mis-ID maps, the biases 
caused by the 'hard'-7r selection of the tagged pion, and 
the application of maps determined in the p region (0.6- 
0.9 GeV/c 2 ) to other mass ranges. These effects are stud- 
ied with simulated 7T7r(7)7isR signal samples, by com- 
paring the mass spectra of produced events when ir- 
ID is either applied or not. The former spectrum is 
obtained by applying the PID process, then correcting 
the '7T7r' spectrum by the 7T-ID efficiency determined as 
in data. The latter spectrum is obtained by not ap- 
plying any PID requirement. As expected, the agree- 
ment is excellent in the p region with a variation of at 
most 2 permil, while some bias is observed in the lower 
and higher mass regions: 1% for to^ < 0.4GeV/c 2 , 
6 x 10~ 3 for 0.4 < m W7r < 0.6GeV/c 2 , 4 x 10~ 3 for 
0.9 < to™ < 1.2GeV/c 2 , and 1% for m OT7r > 1.2GeV/c 2 . 
The full bias determined in simulation is taken as a sys- 
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FIG. 8: The data/MC event correction Cpid for 7tid efficiency 
as a function of m^^. 



tematic uncertainty. The global PID test in data de- 
scribed below supports this estimate. 

In this analysis a good control of the 7r — ¥ '/i' mis- 
ID probability is crucial to the determination of the 
A*A t (7)7iSR cross section in the p region. At track level 
this probability in data is found to vary with momentum 
between 4.5% at 5GeV/c and 7% at 1 GeV/c, and the 
data/MC correction is determined to be 0.95 ± 0.05 and 
1.23 ± 0.03, respectively. 



6. Global PID test with data 

Since the PID classes form an exclusive and complete 
set, every event in the full sample before PID, N xx , is 
assigned to a N%f category ('i','j'= p, tt, K, e, 0). The 
observed N<if spectra and the measured PID efficiencies 
are used in global consistency checks over the full N xx 
data sample. 

The N xx sample is actually composed of Nu pairs of 
particles of identical true types i — p,%,K, with small 
background contributions from other processes that are 
taken into account as follows. The contribution to N xx 
from electrons stems from ecy and 77 followed by a pair 
conversion. It occurs mainly in the l 7re' and 'irir' topolo- 
gies, while being negligible in 'ee' due to the strong re- 
jection of electrons at the preselection and track def- 
inition levels. The small electron component of the 
'7T7r' sample is subtracted out, after proper normalization 
(Sect. IVIB 3[) . In the PID process, protons are mainly 
identified as l ir' and in the global test below, their very 
small contribution is included in N^. Because the N xx 
sample is selected with a tight x 2 ^ < 15 requirement ap- 
plied to the 1C fit, multihadronic background is reduced 
to a negligible level. Contributions from events with two 
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tracks of different true types, from r-pair decays for in- 
stance, are found to be negligible. 

Each observed HV spectrum with 'diagonal' ID, i.e., 
l mr\ l pp\ l KK\ receives contributions from the true 
(ii) channel degraded by the ej_>v efhciencies and from 
the two other channels through ej-}%> mis-ID. The spec- 



tra of produced events in each channel are thus obtained 
by solving a system of three linear equations. In each 
mass bin (computed with the 7T7t mass hypothesis) of the 
spectra for identified pairs of type H\ N>a>, the following 
equations 



are solved for the produced numbers of particle pairs of 
each type, N^^, N^, and N^k- In Eqs. ©, the quan- 
tities Ejj-t<u> represent the product of the ID-efficiencies 
Ej-y'i' and possibly correlation factors that have been es- 
tablished in each PID study. 

From the inferred spectra dNu /dm^n of particle pairs 
of true type i, any Hj' spectrum, dN'if ple ^/dm n „, is de- 
rived, using the measured efficiencies and mis-ID prob- 
abilities, and is compared to the directly observed Hf 
distribution. A relative difference is computed, normal- 
ized to the spectrum dN xx / dm^n of the full sample before 
PID assignment: 

.data _ dN Hrpl . cd /dm^ - dMjf/dm^ 
13 ~ dN xx /dm^ ' 1 j 

All differences <5^ ata are within a few permil. 

The dNxx/dm^jr spectrum is compared to the full in- 
ferred one, dN xx pre d/ dm n7r , obtained by summing the 
N n7r , 7V MM and Nkk components (and the small ee back- 
ground). Figure El which shows their relative difference, 
contains all the information available in data on the va- 
lidity of the ID corrections applied to the different Hi' 
spectra. The band in Fig. [9] represents the limits given 
by the quadratic sum of the estimated systematic uncer- 
tainties on the pp, 7T7T, and KK ID efficiencies. Within 
the statistical uncertainties of the data sample, all devi- 
ations are consistent with the band, thus validating the 
estimates of the systematic errors. 



7. Hard pion specific efficiency 

The hard pion identification is required for one of 
the two tracks when computing the 7r + 7r~ (7)71311 cross 
section off the p peak (Sect. IPX Bl) . The '71V efficiency 
and misidentifications determined in simulation are con- 
trolled in data, and data-MC discrepancies are corrected 
for. The efficiency correction of the '-kit'— > '77717,,' iden- 
tification can only be determined in the central p mass 
region where backgrounds are small in the '7T7r' sample. 
Remaining pp backgrounds are subtracted from the '7T7r' 



and '7r7i7i' samples using the measured e^^' misidenti- 
fication probability and the likelihood selector P^/^ de- 
fined in Section llVD 21 respectively. 

Compared to the standard '7T7r' definition, the event 
ID efficiency is reduced by a factor of 0.825 in data and 
0.870 in simulation. The ratio of the efficiencies in data 
and simulation is shown in Fig. 1101 exhibiting no sig- 
nificant mass dependence between 0.4 and lGeV/c 2 . A 
decrease is observed above 1 GeV/c 2 , which is ascribed to 
an imperfect representation of the large background in 
this region for the '7T7r' sample. A linear fit is performed 
for 0.6 < m n7r < 0.9GeV/c 2 and extrapolated outside 
with propagation of errors. 

8. Separation of the different channels using particle 
identification 

As shown above for the global PID test on data, the 
PID efficiencies and misidentification probabilities mea- 
sured on pure data samples allow one to reliably sepa- 
rate the different two-body ISR channels composing the 
full physics sample before PID assignment. Equivalently 
stated, by solving Eqs. (JSj) on the total physical sam- 
ple, one obtains the produced spectrum dN^/dm^ free 
of misidentified nn{ r y) r fisR and if if (7)71311 backgrounds. 
Likewise, the produced muon spectrum dN^^/dm^^ is 
obtained free of hadronic backgrounds by solving Eqs. © 
per bin of . 

The measured '7r7r' and 'KK' spectra contain con- 
tributions from multihadronic background, from higher- 
multiplicity ISR or qq processes, which are subtracted af- 
ter solving Eqs. ©. There is a contribution from the ISR 
pjry process, which appears dominantly in the '7T7t' spec- 
trum. In this procedure it is treated like pion pairs and 
is removed later. Mistreatment of multihadron events 
where the final state involves a Ktt pair is also consid- 
ered, although they are reduced by PID. 

The procedure of solving Eqs. (JS) to separate the two- 
body ISR channels is applied to isolate the ////(7)tisr 
channel over the full mass range, and the 7T7t(7)7isr chan- 
nel in the p peak region. As the method relies on the com- 
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FIG. 9: (color online). Global PID test on data (see text), 
(top): The m-mr spectrum of all xx^j events (no PID applied, 
data points) compared to dN xx pTed/dm^n (histogram), (mid- 
dle): The different components N xx of the histogram in the 
top plot, with xx labels indicated, and their sum (top his- 
togram with dots), (bottom): The relative difference between 
the two spectra in the top plot: predicted/measured -1. The 
independently-estimated systematic uncertainty is shown by 
the band. 
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FIG. 10: The ratio of the efficiencies for l nnh' identification 
of 7T7t(7)7isr events for data and MC, the line is the linear fit 
result in the p region. 



pleteness of PID class assignment of Table HI it does not 
apply when the '7T7jV selection is required off the p peak 
in the 7T7r(7)7is R channel. In the latter mass regions, 
the A t A t (7)7isR and ^^(t)7isr background subtraction 
is achieved differently as described in Section IVT1 



V. EFFICIENCY STUDIES (II) 

Although the ISR and FSR kinematic fits described 
in Section IIII Bl take into account potential additional 
photons, inadequate description of NLO radiation by the 
simulation might induce incorrect estimates of the 2D-% 2 
efficiency by MC. Comparisons of additional radiation in 
data and MC are performed and data-to-MC corrections 
of efficiencies are applied, as detailed below. 



A. Additional radiation 

1. Additional small-angle ISR 

Additional radiation by the incoming electrons and 
positrons is evidenced by the tail seen along the Xfsr ver " 
tical axis in Fig.[2]for pions and Fig.[3]for muons. In order 
to study the Xisr distributions in data and simulation, 
events are selected above the diagonal (m(x| SR + 1) > 
m (XisR + 1)) an d the restriction E* add ISR > 0.2 GeV 
is applied to ensure a significant level of extra radiation. 
The quantity E* add ISR is obtained from the 2C ISR fit 
described in Section ITlIBl 

The Xisr distributions for the selected events are 
shown in Fig. [IT] As expected for an ISR effect, the 
situation is found to be identical in the muon and pion 
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FIG. 11: Distributions of ln(xf S R + x ) for wMliSR (left) and 7r7r(7)7i SR (right) events in the 0.2-1 GeV/c 2 (0.5-1 GeV/c 2 ) 
mass region for muons (pions), selected with ln(xpsR + 1) > m (xiSR + 1) an d E* a( jd.iSR > 0.2 GeV for data (black points with 
errors) and MC (blue histogram). The distributions are normalized to the data luminosity. 




FIG. 12: Distributions of E* a dd.iSR f° r MA t (7)7iSR (left) and 7T7r(7)7isR (right) events in the 0.2-1 GeV/c 2 (0.5-1 GeV/c 2 ) mass 
region for muons (pions) , selected with ln(xlsR + 1) > hi(xiSR + 1) and El a dd.iSR > 0.2 GeV for data (black points with errors) 
and MC (blue histogram). The distributions are normalized to the data luminosity. 



channels. The data distributions are wider than the MC 
ones. This is a consequence of the ISR fit hypothesis 
that additional ISR photons are collinear to the beams 
as assumed in the AfkQed simulation, as opposed to the 
angular distribution of additional ISR in data. 

The corresponding distributions of E* add ISR , normal- 



ized to the data luminosity, are given in Fig. [T^] for muon 
and pion channels. The MC spectra stop at 2.3 GeV as 
a result of the mx 7 IS r > 8 GeV/c 2 requirement used in 
AfkQed at generation, while the data distributions ex- 
tend to the kinematic limit. Below 2.3 GeV the data and 
MC distributions agree well in shape, but MC is a little 
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higher than data. This is expected since MC includes all 
additional ISR photons, while events in data where the 
additional ISR photon is at large angle have good Xfsr 
and therefore are not present in the sample considered in 
this section. 

The results found for additional ISR in 7T7r(7)7isR and 
fj,fj,("f)"fiSR channels are in agreement, as expected from 
the factorization of additional ISR. The lack of angu- 
lar distribution in AfkQed is studied at the 4-vector 
level using Phokhara, and acceptance corrections esti- 
mated, but its effect essentially cancels in the irir/n/j, ra- 

tio fsect. ivrrrx4i) . 



2. Additional FSR and large-angle ISR 



FSR photons is in good agreement with data and ade- 
quate for our precision goal, since the uncertainty on the 
FSR rate represents about 8 x 10~ 4 of the total /x/x(7)7isr 
sample. 

Figure[T5l fright) shows the same distributions for pions 
after subtraction of the residual 7r + 7r _ 7r°7 background. 
From the comparison of data and MC events in the 
FSR region defined by ln(x FSR + x ) < ln (xfsR + 1), 
E 1 add. fsr > 0.2 GeV, and 9 ni2 < 10°, some excess is 
observed in 7T7t(7)7isr data with respect to the PHO- 
TOS expectation. In the 0.5-1.0 GcV/c 2 mass range, the 
excess is (21 ± 5)%, taking into account subtraction of 
background and the large-angle ISR contribution. This 
difference is accounted for in the determination of the x 2 
selection efficiency, as discussed below. 



Similarly, we select a sample of events with an extra 
photon in the detector acceptance and compares addi- 
tional FSR in data and in the simulation where it is gener- 
ated using PHOTOS. Events are selected with ln(xpsR + 
1) < ln(xisR + 1) an< ^ the fitted energy of the additional 
large-angle photon is restricted to E 1 add. fsr > 0.2 GeV 
in the laboratory frame. The request for a large-energy 
additional FSR photon effectively restricts the 2D-x 2 
plane to a region hi(x 2 SR + 1) > 3. 

The corresponding Xfsr. distributions for /1/X771SR 
data and MC shown in Fig. [T3] (left) are in reasonable 
agreement. Distributions for the selected 7T7T77isr events 
show a similar agreement, Fig. Q2] (right), except in the 
large x 2 tail where there are contributions from back- 
ground in the data. A contribution from secondary in- 
teractions is seen both in data and simulation. 

The components of large-angle ISR and FSR are iden- 
tified looking at the extra photon angular distribution 
with respect to the outgoing charged particles. Figure [HI 
(left) shows the distribution of # M72 , the smallest of the 
two angles with either muon for m M(1 < 1 GeV/ c 2 . A clear 
peak is observed at small photon-muon angle, thus indi- 
cating a true FSR signal in data in agreement with the 
simulation. Evidence for large-angle ISR is also seen in 
data, at variance with AfkQed. This major discrepancy 
is expected, as additional ISR in AfkQed is constrained 
to be collinear with the beams. The same situation is 
observed for pions. After subtraction of a residual back- 
ground from the ISR tt + tt~tt°j process, the FSR peaks 
in data and MC are in fair agreement (Fig. Q3] (right)). 

The photon energy distributions for the /X/Z771SR data 
and MC samples are given in Fig. [TS] (left) for the sub- 
samples satisfying the requirement /i72 < 20°. The 
main component of this sub-sample is from FSR but a 
contribution from large- angle ISR is estimated from the 
(9^ 72 distribution and taken with a 25% systematic un- 
certainty. Absolute rates in data and MC are compared, 
showing a good agreement up to E-y add.FSR ~ 2 GeV, and 
a small excess in data in the tail above. After correction 
for the remaining ISR contribution below 20°, the ratio 
data/MC for the fraction of additional FSR amounts to 
0.96 ± 0.06. Therefore the use of PHOTOS to generate 



B. x 2 selection efficiency 

1. Overview for 7T7t(7)7isr and /i/x(7)7isR 

After the hadronic channels are removed using the PID 
(Sect. IIVD 8|) . the remaining background contributions 
to the jU/i(7)71SR sample are essentially from Trf/y) events 
with two muons in the final state, or one muon and a pion 
misidentified as 'muon'. This contribution is small and 
well simulated, so it can be handled easily even in the 
background (BG) region of the 2D-% 2 plane fSect. lVTA")) . 
Thus the determination of the loose \ 2 selection efficiency 
is straightforward in the muon channel. In contrast, for 
the pion channel, it is not possible to directly measure 
the efficiency of the 2D-% 2 selection in data because of an 
overwhelming background in the rejected region (control- 
ling the loose selection), even in the intermediate region 
(controlling the tight selection). 

The rejected signal events with large % 2 have several 
sources: (i) bad input to the kinematic fits, mostly from 
the direction of the ISR photon, (ii) tails of the % 2 distri- 
butions of events with additional ISR or FSR, (iii) more 
than one additional photon (mostly ISR), and (iv) sec- 
ondary interactions. Except for the last type that is spe- 
cific to pions, the other sources are common to pions and 
muons. A very small difference is also expected for the 
tail of the FSR- fit x 2 , as the FSR level is slightly differ- 
ent for pions and muons. However the level of additional 
FSR is measured in data and MC and the loss of events 
due to FSR can be controlled. 

The strategy is hence to rely on the x 2 selection stud- 
ies performed on muon data to account for the common 
losses and to further investigate the losses specific to pi- 
ons. Therefore the x 2 selection efficiency in data is de- 
rived from the following expression: 

7T7r(7)7isR data _ «t(7)7iSR data X £ 7r // i Qn\ 
X 2 X 2 X 2 ' ' 

where the fj,/ir correction term fe^2 M accounts for two 
effects: i) the difference in additional FSR between pi- 
ons and muons, and ii) pion interactions. The contribu- 
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FIG. 13: Distributions of 1ii(xpsr + 1) f° r A t A i 77iSR (left) and 7T7r77isR (right) events, selected with ln(xlsR + 1) < m (XisR + 1) 
and E 1 add.FSR > 0.2 GeV for data (black points with errors) and MC (blue histogram). The mass regions are chosen to be 
similar for muons (0.2-1 GeV/c 2 ) and pions (0.5-1 GeV/c 2 ); the MC distribution is normalized to the number of events in data 
for muons and with ln(xlsR + 1) < 2 for pions. 



tions of the two components are separated, measured in 
the simulation, and corrected for data/MC discrepancies, 
according to the procedures outlined below. 

2. Determination of the \ 2 selection efficiency for 

The efficiency of the 2D-% 2 selection is measured by the 
rate of hh^jisr events in the rejected region. The spec- 
trum in the signal region is obtained in data by solving 
Eqs. © in each mass bin. The same procedure applied in 
the BG region yields directly the mass spectrum of pro- 
duced muon events rejected by the loose % 2 selection. A 
small contribution from tt must be explicitly subtracted 
using the simulation. The efficiency is directly deduced 
from the ratio of the two spectra. 

Figure [TH] gives the measured x 2 efficiency as a func- 
tion of 171^. It is lower in data than the prediction from 
the simulation by 1.2%. Most of the discrepancy arises 
from the absence of large-angle ISR in AfkQed, which is 
present in data and generates some loss when the large- 
X 2 tails are removed, as expected from Fig. 111! Evi- 
dence for this loss appears at threshold where FSR 
vanishes, while in simulation the efficiency deviates from 
one by only 5 x 1CP 4 . The efficiency decrease with 
is due to the loss of large-% 2 FSR events and reflects 
the increasing rate of additional FSR photons with large 
energy (J5 7 add.FSR > 0.2 GeV). The same behavior is 
observed in data and simulation, consistent with the fact 
that additional FSR is well described in the simulation. 



The systematic uncertainty on the determination of 
the x 2 efficiency comes exclusively from the estimate of 
the background in the muon channel, due to the uncer- 
tainty on the tt cross section and on the PID £j->.y ef- 
ficiencies. These uncertainties amount to a few xlO -4 
and are incorporated in the point-to-point errors. 



3. Effect of additional FSR for pions 

The first component in Se^i^ comes from the difference 
in additional FSR between pions and muons. The differ- 
ence in the FSR rate due to the tt — /i mass difference 
is observed in simulation at the expected level (~ 25%) 
and the expected loss of efficiency due to additional FSR 
is consequently lower by that amount in ^77(7)71311 with 
respect to w^Mtish.- 

The contribution of additional FSR to the \ 2 ineffi- 
ciency is the product of the fraction of FSR events lost 
by the x 2 selection times the rate of FSR events in the full 
pion data sample. To estimate this loss, one relies on the 
Xfsr distributions of the events with additional FSR as 
shown in Fig. 1131 The shapes are similar in the muon and 
pion channels up to ln(x|g R + 1) = 2.5 (which coincides 
with the lower edge of the 'BG' region defined in Fig. 
at which value the pion interaction component turns on. 
The shapes in data agree well with the MC shapes, over 
the entire distribution in /i/x(7)7isR, and, in 7T7r(7)7isR, 
until the distribution in data is affected by background 
in addition to interactions. The fraction of FSR events 
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FIG. 14: The additional 'FSR' photon angular distribution with respect to the closer outgoing muon for the hhjjis'r (left) and 
background-subtracted 7r7T77i S R (right) events with In (xfsr + 1) < In (Xisr + 1)j E 1 a dd.FSR > 0.2 GeV, In (xfsr + 1) < 2.5, 
and in the mass intervals 0.2 < m MM < 1 GeV/c 2 , 0.5 < m W7r < 1 GeV/c 2 (data: black points with errors, MC: blue histogram). 
The mass regions are chosen to be similar for the muon and pion samples and the MC is normalized to data luminosity. 



lost in the pion channel is consequently estimated from 
the fraction of FSR events that fall beyond the Xfsr se- 
lection boundary in the muon data ((35 ± 5)%). This 
fraction of lost FSR events is normalized to the rate of 
additional FSR events observed in the full pion MC sam- 
ple in the ln(xf SR + 1) < 2.5 region (0.64%). The loss is 
further corrected to account for the observed (21 ± 5)% 
data/MC difference in FSR rates (Sect. IVX21) . 

The resulting data/MC correction to the loose \ 2 effi- 
ciency due to the pion-muon FSR difference is estimated 
to be (0.6 ± 0.2) x 10~ 3 . This correction is slightly over- 
estimated as a fraction of the additional FSR events 
is in the signal region of the loose x 2 criterion. The 
correction to the tight \ 2 selection efficiency is larger, 
(1.9 ± 0.8) x 10~ 3 , since all the FSR events with photon 
energy larger than about 0.2 GeV are lost with the tight 
criterion (ln(x 2 sR + -*-) < 3). 



4- Effect of pion interactions 

The effects of secondary interactions are mostly seen in 
the tracking efficiency because of the tight requirements 
imposed on the track pointing to the interaction region. 
The small residual effect in the 2D-x 2 selection efficiency 
is estimated using the simulation, essentially by compar- 
ing the behavior of muon and pion events, and corrected 
for data-MC difference in interaction rates. 

It is found in simulation that the difference of 2D- 
X 2 selection efficiencies, Se^i^, between 7T7r(7)7isR and 



MA*(7)Tisr is about -1.2 x 10~ 3 at 0.75GeV/c 2 . As we 
know that the loss of additional FSR events is smaller 
for pions, since the FSR rate is lower, the smaller effi- 
ciency in 7T7t(7)7isr is ascribed to pion interactions. The 
X 2 selection efficiency loss from secondary interactions 
estimated this way in simulation is (2.8 ± 0.2) x 10~ 3 for 
the loose criteria, flat with mass, and (1.4 ± 0.1) x 10~ 2 
for the tight criteria, with some ±20% relative variation 
with mass. 

Two methods are considered to isolate interacting 
events in both data and MC, and data/MC corrections 
to the x 2 efficiencies obtained above in simulation are 
estimated from the respective rates of pion interactions. 
The corresponding events populate the diagonal region 
in the 2D-x 2 plane, extending through the x 2 selection 
boundary, and therefore affect the x 2 selection efficiency. 

In the first method, interactions are tagged by the pres- 
ence of 'bad' tracks (i.e., tracks not satisfying the track 
requirements of the ISR two-body analysis) in addition 
to the two good tracks of the selected events, provided a 
secondary vertex can be found between a bad track and 
one of the two good tracks. Because of the strict require- 
ments on good tracks, most tagged interactions occur in 
the beam pipe, with further contributions from the first 
SVT layers. The data/MC ratio of interacting events es- 
timated with this method is found to be 1.44±0.10 in the 
intermediate x 2 region, and 1.43±0.13 in the background 
region. However the efficiency of this procedure to tag 
interacting pions is rather low since it keeps about 10% 
of the events with secondary interactions in the signal 
(loose) region and 25% in the background region. 
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FIG. 15: The additional FSR photon energy distributions for the /^77isr (left) and background-subtracted 7T7r77isR (right) 
events with In (xIsr + 1) < m (XiSR + l)i E-y add.FSR > 0.2 GeV, 1h(xfsr + 1) < 2.5, in the mass intervals 0.2 < m MM < 
lGeV/c 2 , 0.5 < m™ < lGeV/c 2 , and # M72 < 20°, #7T7 2 < 10° (data: black points with errors, MC: blue histogram). The mass 
regions are chosen to be similar for the muon and pion samples and the MC is normalized to data luminosity. 



The second method tags a much larger fraction of in- 
teracting events. The quantity doca™* is defined to be 
the largest of the doca xy for the two tracks in the event, 
each limited by the requirement doca xy < 0.5 cm used 
in the good track definition (Sect. The sensitivity 
of this variable to secondary interactions can be appre- 
ciated in Fig. [T71 showing a striking difference in the 
tail of the distributions for pions and muons. The same 
behavior is also observed between pions in events satisfy- 
ing the tight or loose \ 2 criteria. The selection of events 
with doca™y X > 0.1 cm retains about 50% of interactions, 
with a background from non-interacting events that is es- 
timated from the muon distribution. Again it is found 
that the level of secondary interactions is underestimated 
in the simulation, with a data/MC ratio of 1.52 ± 0.03 
in the intermediate \ 2 region. A reliable determination 
cannot be achieved in the background region with this 
method, because of multihadronic background. 

The second determination is more accurate and dom- 
inates the average data/MC ratio of interaction rates, 
which is 1.51 ± 0.03. 



5. Summary of corrections to the \ 2 selection efficiency for 

pions 

Adding the two components of the \ 2 selection ineffi- 
ciency in 7T7t(7)7isr that are not common to the muon 
channel, the extra data/MC correction to apply to <5e^2 P 
for the loose \ 2 selection amounts to (1.4 ± 0.1) x 10~ 3 



for secondary interactions and (0.6 ±0.2) x 10 3 for FSR. 
The total correction is (2.0 ± 0.3) x 10~ 3 . 

For the tight \ 2 selection, both corrections are larger: 
(7.1 ± 0.4) x 10" 3 for interactions and (1.9 ± 0.8) x 10" 3 
for FSR. The total data/MC correction on the tight \ 2 
selection efficiency amounts to (9.0 ± 0.9) x 10~ 3 . 



VI. BACKGROUNDS 
A. Backgrounds in the /i^(7)7isr channel 

Contributions from hadrons to the '/i/i' sample are re- 
moved by PID through solving Eqs. ([8]), as described in 
Section IIV D 8[ but one must still consider background 
from processes producing real muons. 

ISR-produced J/ip followed by decay to /i/i is not a 
background to the complete /x/x(7)7isr process, but is 
a background to the purely QED reaction used for the 
determination of the ISR luminosity. The '0(25') con- 
tributes as a background through its decays to J/ip, ei- 
ther following the ir Tr°J/ip transition or radiative decays 
through charmonium states. Both direct and indirect 
J/ijj production is observed. These contributions are re- 
moved by excluding events where the measured invariant 
mass is in the 3.0 — 3.2 GeV/c 2 window 4 . 



4 In the final mass spectrum, this rejection does not however 
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FIG. 16: The \ 2 efficiency (top) for /i/u(7)7isR data (after 
background subtraction), MC (AfkQed), and the ratio C x 2 
of data to MC (bottom), as a function of m MM . The gap 
at 3.0 — 3.2 GeV/c 2 corresponds to the excluded J ftp window 

(Sect. East . 



Another hidden background from J/ip comes from the 
radiative decay J/ip — ¥ /x + ^ _ 7, which is indeed barely 
observed in the ///iry mass spectrum in data. Its contri- 
bution to the n/j, mass spectrum between 2 and 3 GeV/c 2 
is of order 10~ 3 and neglected. 

The process e + e~ — > t + t~(7) can contribute to the 
'/z/k' sample through r — > pv ll y T and mis-identified r — > 
7w T decays. The contribution is estimated by MC and 
found to be negligibly small, except at masses above 



produce a sharp hole as m MM is determined after the (J-fJ-i^fiSR 
kinematic fits described in Sectior flll Bl 




0.1 0.2 0.3 0.4 0.5 
docaxy(cm) 



FIG. 17: (top): The distribution of the largest of the two 
transverse distances of closest approach to the interaction 
point doca^y* for pions and muons in data, for the interme- 
diate \ 2 region, (bottom): Same distributions in simulation. 

2 GeV/c 2 where it reaches a fraction of 10 ~ 3 . It is sub- 
tracted using the simulation. 



B. Backgrounds in the 7T7t(7)7isr channel 

1. Background from /!/i(7)7isr. 

Separation of each component of the two-prong ISR 
sample is achieved through solving Eqs. © in each m n7r 
mass bin. This procedure yields the produced spectrum 
dN n ir / d-m^TT and the background contributions to the ob- 
served '7r7r'-identified spectrum. The muon background 
level is less than 4 x 10~ 3 at the p peak but increases 
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rapidly away from the resonance and reaches a few per- 
cent at the p tail boundaries (Table HIT)) . 

As the '7T7tV selection applied in mass ranges away 
from the p precludes using the above procedure, the re- 
duced /u/h(7)tisr background contribution is determined 
directly from the L nirh sample using the track identified 
as a l 7r' with the standard pion identification. The fit 
of the distribution, for that track, of the n/fi likelihood 
estimator (introduced in Section llVD 2 1) yields the 
respective true muon and pion components of the 'tt'- 
identified tracks, hence of the '^717^ sample. Fits are per- 
formed in 0.5 GeV/c 2 -wide mass bins and, except for the 
0.5-1 GeV/c 2 interval (p peak region), the pp component 
can be well determined. Above 3 GeV/c 2 the muon con- 
tribution becomes dominant, despite the '7r7i"h,' ID, and 
the pion signal is lost above 4GeV/c 2 . The results of 
the fits are summarized in Fig. [TSl which shows, for each 
mass interval, the pp — > '7T7tV fraction in data relative 
to the prediction from the simulated muon sample after 
luminosity scaling. A second-order polynomial fit to all 
points allows one to smoothly interpolate between the 
low and high mass regions. The band indicates the error 
envelope of the fit. 
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FIG. 18: The rate of pp events misidentified as '7T7iV in data 
relative to the MC expectation, measured in mass intervals 
below and above the p region . The curve with the error 
band is a second-order polynomial fit to the data points, used 
to interpolate through the p mass region (0.6 — 0.9GeV/c 2 ). 



2. Background from K K {7)71311 

When the standard 7T-ID identification is applied to 
both tracks in the '7T7r' sample, the kaon background is 
implicitly subtracted through solving Eqs. © and stays 
below the permil level in the p peak region (Table IIIII) . 



This background is essentially insensitive to the fur- 
ther selections applied to the '7T7r^' sample. Since the 
.?C-R"(7)7isr events are dominated by the narrow <f> res- 
onance, one can use this feature to determine the KK 
component directly in data. The procedure is illustrated 
in Fig.1191 which shows the ttikk mass distribution of the 
'7T7r^' sample when the K mass is assigned to both tracks. 
A <f> signal is fitted with the signal line shape taken from 
simulation and a linear term to describe the dominant 
true 7T7r component. The <fi signal yield provides the nor- 
malization of the remaining K K (7)7isr background con- 
tribution to be subtracted from the 'iri^h data. The wide 
and distorted shape of the <f> peak reflection in the m nv 
mass spectrum is taken from simulation. 
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FIG. 19: The ttikk mass distribution in the '^^^7' data sam- 
ple. The solid line represents the result of the fit (see text). 



3. Background from ee"/ events 

Radiative Bhabha events are very strongly suppressed 
in the event selection because of the track definition that 
contains a veto on electrons. Remaining events of this 
type are from distribution tails and various pathologies. 
Because of this large selection bias there are very few 
events actually identified as 'ee'7 in the identification 
process. Radiative Bhabha background appears in the 
'en - ' and '7T7r' ID topologies. This background is iden- 
tifiable near threshold and at high masses. Its small 
contribution cannot be detected in the p region and its 
mass shape and magnitude are estimated from a down- 
scaled sample of radiative Bhabha events normalized near 
threshold. 

The radiative Bhabha background normalization is 
achieved using the angular distribution in the 7T7t center- 
of-mass system, assuming the pion mass for the parti- 
cles. The angle 6* of the ir + is measured relative to the 
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ISR photon direction in that frame. In the mass range 
0.28 < m wv < 0.32GeV/c 2 just above threshold, the 
remaining eej background contribution, still noticeable 
with the 'iTTTh identification, is obtained by fitting the 
| cos 9*\ distribution. Backgrounds from pp^jisn and 
KK^jisr are subtracted before fitting, with shapes 
taken from simulation with correction from the data, 
and normalized to data luminosity. The | cos 9* | distri- 
bution is fitted with two components: mr^^/iSR-i with 
the shape taken from the simulation, and ee7 background 
with the shape obtained from the downscaled radiative 
Bhabha sample. The latter contribution has a charac- 
teristic sharp peak near one with a long tail while the 
7T7r(7)7isR signal behaves as sin 2 0*. The fit shown in 
Fig. [20] on the '71717,' sample at threshold provides the nor- 
malization factor to be applied to the radiative Bhabha 
sample to describe the eej background. A similar fit per- 
formed at threshold on the l irir' sample yields a consistent 
normalization factor. 




FIG. 20: (color online). The | cos 6 1 * | distribution of l TVKh 
data in the 0.28-0.32 GeV/c 2 m^^ range, fitted (black curve) 
to two free components: 71^(7)773^ from MC (blue dashed 
line) and ee7 background from downscaled radiative Bhabha 
events (red dotted line). The small pp{~f)^isR contribution 
is subtracted out. 



for m wn < 0.32GeV/c 2 and less than 0.05% for larger 
masses for the 'irirh selection, and 0.63(0.27)% at m n7r = 
0.525(0.975) GeV/c 2 for the 'tttt' selection. At the p peak 
the fraction drops to only 0.03%. 

A systematic uncertainty of 100% is assigned to the 
ee normalization factor determined at threshold, which 
is applied up to 3 GeV/c 2 . This precision is adequate in 
the p region as well as in the resonance tails. 



4- Conversions and rejection of displaced vertices 

Purely electromagnetic processes may induce back- 
grounds when one of the final-state particles interacts 
with the detector material, allowing the selection crite- 
ria to be satisfied. This is the case at threshold from 
the e + e~ — > 77 process followed by a photon conversion, 
and at large masses from Bhabha scattering where one 
of final electrons (positrons) undergoes bremsstrahlung 
in the beam pipe. In either case, one or both of the 
detected tracks can be misidentified as pions. However, 
as they do not originate from the interaction point, this 
contamination is reduced by requiring the distance in the 
transverse plane V xy of the vertex of the two tracks and 
the average interaction point to be small. 

The background from conversions is expected to yield 
a rather wide V xy distribution, while prompt particles 
(from ee7 and ^71(7)7/5^) produce a peak at zero. The 
requirement V xy < 0.5 cm is applied in the p tails, which 
are the mass regions affected by the background from 
conversions. The conversions are reduced to a negligible 
level (< 5 x 10~ 4 ) by the selection. The efficiency of the 
Vxy requirement for ^77(7)7/5^ events is controlled by 
the two-track vertexing and pion secondary interactions. 
The former effect is studied in data and simulation using 
the ppi^jiSR sample. The effect of the pion secondary 
interactions is studied in the p region, both in data and 
MC. The overall correction to the MC efficiency is (1.1 ± 
0.1)% at 0.4GeV/c 2 and smaller for larger masses. 



PP7isr process 



The mass dependence of the eej background is checked 
using the sample of events identified as 'en'. The latter is 
a rather pure eej sample outside the p region, and is more 
representative of events mis-identified as 'mr^ than the 
genuine radiative Bhabha downscaled sample. The ratio 
of the mass spectra of 'e7r' events to the radiative Bhabha 
sample is found to be constant within uncertainties, for 
masses away from the p peak, i.e., just above threshold 
and in the 1.5-3 GeV/c 2 range. 

The ee7 background is only noticeable at thresh- 
old and near the edges of the p central region: 6% 



Proton ID is not considered in the particle identifi- 
cation process, since the process pjryisR contributes at 
a very small level. With the chosen ID classes pro- 
tons are classified as pions, and antiprotons sometimes 
as electrons. The cross section for the pfryiSR process has 
been measured by BABAR [22[ and the results are used to 
reweight the MC prediction. The overall contamination 
is taken from the reweighted simulation and subtracted 
statistically. It amounts to less than 0.5% in the p central 
region (Tab. IIIII) and exceeds the percent level at large 
masses only (m vn > 1.1 GeV/c 2 , Tab. HV|) . 
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6. Multihadrons from the qq process 

Hadronic processes, either direct or ISR-produced, in- 
troduce a background in the pion sample that is con- 
siderably reduced by the requirement of only two good 
tracks and the x 2 selection of the kinematic fits. This 
contribution is estimated using simulated samples of the 
e + e _ — > qq process. However, the JETSET prediction for 
qq fragmentation into low-multiplicity final states is not 
necessarily reliable, so the MC rate is normalized using 
data. 

In backgrounds from the qq process, the ISR photon 
candidate actually originates from the decay of an ener- 
getic 7T°. Such a signature is searched for, both in data 
and in MC, by pairing the ISR photon candidate with all 
detected additional photons and the MC normalization 
is obtained from the observed 7r° rates. The pair with 77 
mass closest to the nominal tt° mass is retained. Fits to 
77 mass distributions are performed in data and MC as- 
suming a Gaussian shape for the ir° signal and taking into 
account 'background' from 7T7t(7)7isr events and contri- 
butions from other processes, like e + e~ — > tt(j) events, 
both taken from the simulation. 

The 7T° fits are carried out in wide irir mass bins 
(0.5GeV/c 2 ) between threshold and 3GeV/c 2 , covering 
the practical range for the analysis. Fits are also per- 
formed in background enriched regions to check the sen- 
sitivity to the final-state multiplicity. An example of fits 
is shown in Fig. [5T] JETSET is found to overestimate 
the background contributions by a factor of 1.3, almost 
independent of the 7T7t mass and whether it is determined 
in the signal or background enriched regions. 




nyGeV/c') 

FIG. 21: The 7isr7 mass distribution for data events in a 
background enriched region. The ir° signal is fitted with a 
Gaussian while additional contributions are represented by a 
linear term. The histogram is the 7T7t(7)7isr MC distribution. 
Contribution from rr events has been subtracted. 



We check that the ir° finding efficiency does not depend 
on the final state produced by fragmentation in JET- 
SET. The results are consistent within 5% for ir + ir~ir°, 
7T + 7r~27r°, and the full qq contribution. The final state 
with lowest multiplicity, 7r + 7r~7r°, which is topologically 
and kinematically identical to the signal, is further con- 
trolled on data. As the 7T7t mass distribution for these 
events in JETSET peaks between 1 and 2 GeV/c 2 , well 
above the p peak, we search for a tt° signal in that mass 
range where the e + e~ — > 7T7t(7)7isr contribution is small. 
The weak signal observed in data is consistent with the 
negligible JETSET expectation within a 50% uncertainty. 

After renormalization, the qq background fraction 
is below 1.8% in the loose \ 2 region (Tab. and 
much smaller when the tight \ 2 criterion is applied, ex- 
cept above 1.3GeV/c 2 where it reaches a few percent 
(Tab. EE). 

The MC statistical errors are included in the sub- 
tracted qq background spectra, while the uncertainty in 
the normalization (from 7r° fits), typically 10%, is taken 
as a systematic error. 



7. Multihadronic ISR processes 

The background is estimated using simulated processes 
e + e _ — !• A71SR where X stands for the final states: 
7r + 7T~7r°, 7T + 7r~27r , 27r + 27r~, 27r + 27r~7r°, r/7r + 7r~, and 
KsKl. The dominant contributions are from e + e~ — ¥ 
7r + 7r~7r°7 and e + e~ — > 7r + 7r~27r°7. They sum to about 
10% at the lower edge of the p peak (Tab. IHIj) but are 
strongly reduced by the tight \ 2 criterion. 

An approach similar to qq comparing data and MC is 
followed for the ir + ir~ir° ISR normalization. This process 
is dominated by the production of the to and <j) narrow 
resonances, which are used as calibration signals. A kine- 
matic fit to the 7r + 7r~37 final state is performed using a 
7T° constraint, and the 7r + 7r~7r° mass distribution is fit- 
ted. The ratio of the contributions in data and MC is 
found to be 0.99 ± 0.04. The error on the normalization 
factor is taken as a systematic uncertainty. 

The remaining ISR processes are higher-multiplicity 
7T° hadronic states such as 27r37r°7. These cross sections 
have not yet been measured by BABAR but we estimate 
that the contributions of these channels to the total back- 
ground do not exceed the 10 -3 level. It is estimated from 
MC alone, assuming a normalization uncertainty of 10%. 

8. Background from other processes 

The e + e~ — > T + T~("f) process contributes significantly 
only at 7T7T masses higher than the range of interest for 
this analysis. Although at a very small level, this back- 
ground is subtracted using simulation. Two-photon pro- 
cesses with hard radiation such as e + e~ — > (e + e~)7r + 7r~7 
and the similar reaction with muons have been specifi- 
cally looked for in kinematic regions where they are ex- 
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pected to contribute, but without finding a significant 
effect. 

A summary of backgrounds and related errors are given 
in Tables HI] and [TV] 



9. Overall test of the multihadronic background 

The multihadron background fraction estimated above 
reaches sizeable values near the boundaries of the central 
p region, but with a quite small uncertainty, 4.8 x 10~ 3 
at 0.5GeV/c 2 and 3.0 x 1(T 3 at 1GcV/c 2 , the value at 
the p peak being negligible. In the p tails, the estimated 
systematic errors due to multihadron backgrounds, which 
are strongly reduced by the tight \ 2 criterion, do not 
exceed a few permil. 

We assess both the rate and the mass distribution of 
the multihadron background in data, in the 2D-x 2 re- 
gion where it is the largest, i.e., in the 'sleeve' outlined 
in Fig. [22] We fit the ttit mass distribution in this re- 
gion (Fig. (23]) to background and signal components, with 
shapes taken from MC. The ratio of the fitted background 
to the one estimated in the above sections is found to 
be 0.968 ± 0.037, consistent with unity. This is trans- 
lated into an uncertainty of 4.5 x 10~ 3 at 0.5 GcV/c 2 and 
1.5 x 10~ 3 at 1 GeV/c 2 on the background fraction in the 
full 7T7T sample. These values are below the quoted un- 
certainties, which validates the multihadron background 
estimate and confirms that the mass distribution of the 
background from the simulation is appropriate. 



qqMC 0.5< miIIt <1.0GeV/c 2 




10 

ln (Xl S R+l) 



FIG. 22: (color online). 2D~x 2 distribution of qq MC events 
normalized to the data luminosity for 0.5 < ni^-n < 1 GeV/c 2 . 
The solid broken line indicates the loose \ 2 criterion, while 
the dashed line defines a 'sleeve' in the signal region where 
most of the background is concentrated. 



> 



4000 - 



c 

u 
> 



2000 




m^GeV/cf) 



FIG. 23: The m^^ distribution in the background-rich 'sleeve' 
region. The solid line represents a fit to the data with ttit 
signal and multihadron background components, with their 
shapes both taken from simulation. 



VII. MASS SPECTRA DETERMINATION 

The spectra of pp("j)jiSR and 7^(7)71311 events after 
event selection are obtained as functions of the two-track 
mass (m„) given by the best x% (Xnn) nt - These 
spectra are background-subtracted, and mass-dependent 
corrections for data/MC efficiency differences are applied 
as described above. To account for FSR effects and res- 
olution smearing due to the detector response, unfolding 
is required to obtain the dN/d^fs' spectra as functions 
of the final-state mass including FSR, which are used to 
measure the cross sections through Eq. ([I]). 



A. Unfolding of the mass spectra 



1. The unfolding method 

The unfolding technique used in the present analysis 
is a simplified version of a method developed for more 
complex unfolding problems [23] . The folding probability 
Pij of an event produced in a true (Vs*) bin j to be 
reconstructed in a (m xx ) bin i is computed directly in 
simulation from the transfer matrix Ay (the number of 
events produced in a true bin j that are reconstructed in 
bin i) 5 . Conversely, the matrix of unfolding probabilities 



The matrix of folding probabilities is related to the transfer ma- 
trix Aij by Pij = Aij j X)fe=i ^fej while the matrix of unfolding 



probabilities is 
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PL indicates the probability for an event reconstructed 
in a bin i to originate from the true bin j, and is also 
computed from the transfer matrix. Ay and P[a depend 
on the assumed true spectrum while Py- , which describes 
detector and FSR effects, docs not. The method used to 
unfold the m xx spectra is based on the idea that if the 
MC describes well enough the true spectrum in data and 
if the folding probabilities are well simulated, the matrix 
of unfolding probabilities determined in simulation can 
be applied to data. 

If the first condition is not fulfilled, that is if the data 
spectrum after unfolding differs significantly from the 
true MC spectrum, several steps are considered where 
the transfer matrix is improved by re- weighting the true 
MC, keeping the folding probabilities unchanged. Dif- 
ferences between data and folded ('reconstructed') MC 
spectra are ascribed to differences in the unfolded ('true') 
spectra. At each step of the iterative re-weighting pro- 
cess, the data-MC differences of reconstructed spectra 
are unfolded and added to the true MC spectrum. Such 
iterative procedures can result in a significant bias to the 
final results if statistical fluctuations are mis-interpreted 
as true differences between data and MC distributions. 
The stability of the method is provided in this analysis 
by the use of a regularization function to avoid unfolding 
large fluctuations in the data, due for example to a large 
background subtraction. Details on the method are given 
in Ref. [Uj]. 



2. Procedure 

Unfolding is applied to the reconstructed m xx spec- 
trum, after background subtraction and data/MC cor- 
rections for efficiencies, obtained as described in the pre- 
vious sections. The unfolding procedure handles detector 
resolution and distortion effects, and corrects for FSR. 
Thus the process delivers the 'true' distribution of events 
in the detector acceptance as a function of s/s', and a 
covariance matrix containing the statistical uncertainties 
and their bin-to-bin correlations. The covariance ma- 
trix is obtained from pseudo-experiments (toys), where 
both the spectrum and the transfer matrix are statisti- 
cally fluctuated. 

For the 7T7t(7)7isr analysis, the same energy range 0- 
3 GeV is chosen for data and the MC transfer matrix. The 
spectra obtained under the central p region conditions 
(loose x 2 criterion) and under the p tails conditions (tight 
X 2 criterion) are unfolded separately over the full mass 
range. The unfolded spectra are combined afterwards, 
each being used in its respective mass region. Different 
bin sizes are used: 10MeV/c 2 for the tails and 2MeV/c 2 
for the central part. 



Since the loose condition retains events in the inter- 
mediate 2D-x 2 region, resolutions in data and simulation 
are compared in specific ranges across this region. It is 
noteworthy that, although the large-angle additional-ISR 
events, absent in MC, populate the intermediate region, 
they do not contribute to resolution tails in data as the 
mass in that case is given by the (good) FSR fit. Mass 
spectra shapes are found to be well simulated, but the 
rate of events in the degraded mass resolution regions is 
underestimated by MC. A reweighting of the MC sample 
in the corresponding x 2 regions is applied, thus modify- 
ing the MC transfer matrix. The relative distortion of 
the mass spectrum due to reweighting is less than one 
percent (see Fig. HIT) , hence the systematic uncertainty 
from the imperfect knowledge of the transfer matrix is 
estimated to be a fraction of 10~ 3 . 

The unfolding of the monotonic featureless mass dis- 
tribution in the Hh(j)jisr sample is much less sensitive 
to resolution effects and the transfer of events involves 
mostly FSR. Here a larger mass range (0-6 GeV/ c 2 ) is 
considered in 50MeV/c 2 intervals (120x120 matrix), al- 
though only the first half is needed for luminosity pur- 
poses. 

The initial mass-transfer matrices for PP^Jisr 
events, and 7T7r(7)7isR under the loose and tight condi- 
tions, are shown in Fig. [531 The large diagonal compo- 
nent corresponds to a mass resolution 6 of 3.2MeV/c 2 at 
m n7r = 0.78 GeV/ c 2 (loose x 2 criterion), and 4.6MeV/c 2 
at m W7T = 1.5 GeV/ c 2 (tight x 2 criterion). 

The most significant difference between data and re- 
constructed MC in relative terms occurs in the region 
1.7-2 GeV/c 2 , where the pion form factor is not well sim- 
ulated. Smaller differences, not exceeding the statistical 
errors for 2-MeV bins, are observed in the p lineshape, 
in the tails and in the peak region with the p — uj in- 
terference (Figf55] (top)). These differences are assigned 
to the generated mass distribution in the MC, as resolu- 
tion effects between data and MC are studied separately: 
the resolution broadening in the intermediate % 2 region 
was discussed above and the resolution for the tight x 2 
condition is presented in Section IVII CI The differences 
are corrected for in the iterative way described above, 
but it is observed that already after the first step of the 
re-weighting procedure, they are reduced to a negligi- 
ble level. The residual systematic differences have in- 
deed very little effect on the result of the unfolding. The 
first unfolding result is very close to the initial data (well 
within the statistical error), except in the p — uj inter- 
ference region, as expected since the mass resolution is 
not small compared to the to width (8.5 MeV). Adding 
one iteration in the unfolding does not result in further 
improvement, as shown in Fig. 1251 (bottom). 



6 For illustration here we use an effective mass resolution obtained 
of events produced in a true bin j that are reconstructed in bin by taking the weighted-average of the standard deviations from 

i. a two-Gaussian fit of the resolution function in simulation. 
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FIG. 24: The initial mass-transfer matrix Aij from the simulation giving the number of events generated with a (true) mass 
y/s' in a bin j and reconstructed with a (measured) mass m xx in a bin i: (J,fJ,{ , y)'yisR (top), 7T7t(7)7isr with tight \ 2 criterion 
(middle) and with loose \ 2 criterion (bottom). For the latter case only the relevant range 0.5-1.0 GeV/c 2 is shown. The y/s 1 
dependence comes from QED and a model of the pion form factor used in the AfkQed generator, respectively. 



3. Tests of the unfolding method 

Tests of the unfolding procedure are performed, in the 
fifi and tttt channels, investigating potential systematic 
biases introduced by the method. The test uses toy dis- 
tributions of true spectra and their corresponding recon- 
structed distributions obtained by folding using the nom- 
inal transfer matrix A. The reconstructed toy spectrum 
is then unfolded with a transfer matrix (A) obtained after 
statistically fluctuating A. Finally the unfolded result is 



compared to the true toy spectrum. 

The true toy distribution is constructed from the true 
MC with a bias added. In order to build a test as close 
as possible to the real situation, the bias is taken as the 
difference between data and the normalized initial recon- 
structed MC. Two variations of the test are considered, 
where the reconstructed spectrum is additionally fluctu- 
ated statistically or not. The first situation is closer to 
the real unfolding operation and could reveal spurious ef- 
fects due to the limited statistics in the data (and MC). 
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FIG. 25: (color online), (top): The difference between the 
tt + tv~ mass distributions of data and reconstructed MC before 
unfolding (data— rMC) and after one iteration (data— rMCm), 
for the loose \ 2 selection used in the central p region. The 
data statistical errors (±lcr) are shown for comparison. The 
correction to the initial MC distribution is small, but signif- 
icant in the peak and tail regions, (bottom): The difference 
between the result of the first unfolding (UR1) and the initial 
data for the same loose x 2 criterion. It exceeds the data sta- 
tistical error (band) only in the p — uj interference region. No 
significant improvement is observed between the first (UR1) 
and second (UR2) unfolding results. 



The second test allows one to search for potential sys- 
tematic effects of the method itself. 

It is found that the systematic bias on the yfs' spec- 
trum from the unfolding technique is negligible in the 
/jL/i channel. In the mr channel, it is below the 10~ 3 
level, except in the ranges 0.5-0.6 GeV (1.9 x 10 -3 ) and 
0.9-1.0 GeV (1.2 x 10" 3 ). The latter two values are anti- 
correlated with the rest of the spectrum, hence the sys- 
tematic uncertainty that affects the dispersion integral 
used in the a„ calculation remains smaller than 10 -3 . 



persion integral well below 10 



B. Mass scale calibration 

The absolute 7T7r mass scale depends on the momenta 
and angular measurements and the kinematic fit. Un- 
like at threshold where the mass scale is governed by the 
angular measurements, the uncertainty from the momen- 
tum scale is dominant at the p mass and above. There- 
fore systematic effects are studied using ISR-produced 
J/ip — > fifj, events, which are treated in the same way as 
the di-pion sample. 

The /i/i mass distribution is fitted in the 3.0-3.2 GeV/c 2 
range across the J / ip peak with a linear term for the QED 
background and a signal shape obtained by convoluting 
the sum of the natural J / ip Breit-Wigner and the QED- 
J/ip interference with a Gaussian resolution shape. The 
free parameters are the amplitude of the signal, the J/ip 
mass mji^,, the resolution a m and the two constants of 
the linear background term. Three bins (1-3; 3-5; 5- 
8 GeV/ c) are defined for the two track momenta p m m and 
< Pmax) and fits are performed in 6 boxes in 
the (pmin,Pmax) plane, not distinguishing charges. 

Whereas a m increases for larger momenta as expected, 
the fitted values for mj/^ are consistent for all boxes, 
showing no evidence for a momentum-dependent calibra- 
tion change. Therefore the whole sample is considered, 
the corresponding fit being shown in Fig. 1261 The small 
symmetric excess observed in the tails does not affect 
significantly the central value. The result 



m J/4 , = (3096.30 ± 0.13) MeV/c^ 
a m = (9.38 ± 0.04) MeV/c 2 , 



(11) 
(12) 



is compared with the world-average value 24] , (3096. 92± 
0.01)MeV/c 2 . The difference, (-0.62 ± 0.13) MeV/c 2 , is 
interpreted as a momentum scale shift of (— 2.00±0.04) x 
10~ 4 . 

This momentum-calibration scale factor translates into 
a shift for the p mass of (— 0.16±0.16) MeV/c 2 where the 
full correction is taken as a systematic uncertainty. 



4- Distortion of the mass spectrum due to excess FSR 

The small excess of events with additional FSR in 
data compared to the simulation produces a distortion of 
the mass spectrum not taken into account in the mass- 
transfer matrix. By appropriately reweighting the energy 
distribution of FSR photons by the energy-dependent ex- 
cess fraction one obtains the resulting systematic uncer- 
tainty on the mass distribution. The maximum deviation 
in the p region occurs at 0.5-0.6 GeV/c 2 at the 2 x 10~ 3 
level, while it decreases to —0.8 x 10~ 3 at the p peak 
and —0.5 x 10~ 3 at 1 GeV/c 2 . These values are taken as 
systematic uncertainties on the cross section. Because of 
the anti-correlation occuring below and above the peak, 
this effect produces a systematic uncertainty on the dis- 



C. Mass resolution 

Since detector resolution effects are corrected for 
through the unfolding procedure using the mass transfer 
matrix from simulation, we check that simulation repro- 
duces data in this respect. This problem is not crucial 
since the 7r + 7r~ annihilation cross section is dominated 
by the wide p resonance, except in the p — uj interference 
region because of the small w width (8.5 MeV compared 
to the 7.6 MeV FWHM resolution). 

In the p region the mass resolution is dominated by mo- 
mentum, rather than angular measurements. Thus the 
J/ip study described in Section I VII Bl is again relevant. 
The mass resolution found in data (Eq. (fT2"j)) is slightly 
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FIG. 26: Fit of the 'fj,fj,' mass distribution in the J/ip region 
including the QED- J/ip interference as a momentum calibra- 
tion test. 



better than the result from the simulation of continuum 
w(l)jiSR events in the 3.0-3.2 GeV/c 2 range (no J/ip 
contribution is generated in AfkQed), which is found to 
be (10.0 ± 0.1) MeV/c 2 . 

The contribution of the decay angle measurement to 
the mass resolution is obtained for data and simulation 
from a study of the decays Kg — » ir + ir~ , from a sam- 
ple of ISR- produced <p mesons decaying into KgK®. In 
this case the angular measurement plays the dominant 
role compared to momentum, thus this is a situation 
complementary to the J/ip one. Taking into account 
the smaller contribution from momentum taken from the 
J/ip case, the study yields the average resolutions on the 
decay opening angle of (1.65 ± 0.03) mrad in data and 
(1.59 ± 0.03) mrad in simulation. 

Combining the momentum and angular contribu- 
tions one obtains the full mass resolutions of (3.03 ± 
0.03) MeV/c 2 and (3.20 ± 0.03) MeV/c 2 for data and sim- 
ulation, respectively. This resolution difference results in 
a bias on the measured resonance widths after unfolding 
of the mass spectrum, given by 

Ar p = (+0.016 ± 0.004) MeV, (13) 
Ar u = (+0.27 ±0.07) MeV. (14) 

As for the mass scale calibration the full biases are taken 
conservatively as the corresponding systematic uncer- 
tainties on the measured p and u> widths. 



VIII. RESULTS ON e + e~ /i + M~(7hr3R CROSS 
SECTION AND ISR LUMINOSITY 



Simultaneous measurement 
V + V~ (l)lisn and e+e~ 



of the e + e — > 
r + 7r~ (7)7isr channels 



is a major feature of this analysis. In this section, we 
report the results on the absolute e + e~ — > /i + M~(7)7iSR 
cross section measurement and the comparison to QED. 
The measured mass spectrum is corrected for 

all efficiencies described in the preceding sections and 
unfolded. Further corrections are specific to the absolute 
cross section measurement, which necessitate dedicated 
studies described in this section. We then express 
the results on the e + e~ — > (7)7isr spectrum 

in terms of the effective ISR luminosity used in the 
e + e~ — » 7r + 7r~(7) cross section measurement. 



A. Acceptance effects specific to the pp{'y)'yisR 
analysis 

1. Relevance of these studies 

As stressed in the Introduction, the measurement of 
the e + e~ — > cross section relies on the measure- 

ment of the mr/n/i ratio R exp (y/~s')- A major advantage 
is that the effect from additional ISR essentially cancels 
in the ratio, leaving only second-order effects that are 
studied specifically. 

However, for the QED test, we use the absolute mea- 
surement of the ^A t (7)7isR cross section in order to per- 
form a direct comparison to QED at NLO. This is a strin- 
gent check of the full understanding of all involved sys- 
tematic effects. Dedicated studies are conducted in order 
to assess the importance of NLO effects in the MC gen- 
erator. The QED test, if successful, demonstrates that 
these effects are properly taken into account and that 
their residual impact on the irir/fifi ratio measurement, 
in which they largely cancel, is indeed very small. 



2. Extra radiation in the MC generators 

We use AfkQed as the h/i^^isr event generator 
to produce a large sample (5.3 times the data) of 
fully-simulated events. The radiator function dW/d\fs' 
and vacuum polarization correction a(s')/a(0) entering 
Eq. ([6]) are included in the generator. However, although 
AfkQed describes correctly the lowest-order process, it 
has some shortcomings in the generation of extra ra- 
diation: (i) additional ISR photons are generated with 
the structure function method in the collinear approx- 
imation, with a photon energy cut-off near 2.3 GeV in 
the e + e~ cm. (coming from the requirement rnx~f lSTi > 
8 GeV/c 2 applied at generation) , (ii) generation of addi- 
tional FSR photons follows the PHOTOS algorithm. 

The effects of these limitations are studied with the 
Phokhara 4.0 generator [l7j]. The advantage of Phokhara 
is that it uses the almost-exact QED NLO calculation 
(without ISR-FSR interference). However, contrary to 
AfkQed, it does not include the contribution from two 
FSR photons nor higher-order ISR emission. Both effects 
are expected to be at a very small level. The contribution 



33 



of two FSR photons is suppressed by the smallness of LO 
FSR (about 1% at lGeV and 15% at 3GeV) and NLO 
FSR (< 1% at 1 GeV and 2.7% at 3 GeV) corrections for 
photon energies E 1 add.FSR > 0.2 GeV. Even at 3 GeV 
the expected contribution of 4 x 1CP 3 has a negligible 
effect on the acceptance. The contribution of higher- 
order ISR emission is relevant only if the third photon 
has a significant energy. From the acceptance change 
between Phokhara and AfkQed, and the fraction of NLO 
ISR above photon energies of 1 GeV in the e + e" cm., 
one estimates a maximum acceptance bias of 2 x 10~ 3 at 
threshold and 1CT 3 at 1 GeV. 



3. Fast simulation studies with Phokhara and AfkQed 

Since any number of additional photons are accepted 
at event selection, an imperfect simulation of NLO ISR 
does not affect the event topology selection but alters the 
event acceptance through kinematic effects. 

The main criteria affecting the geometrical acceptance 
are: both muon tracks in the polar angle range 0.4 < 
Ofj, < 2.45 rad, with momenta larger than 1 GeV/c; the 
most energetic photon in the cm. (ISR candidate) with 
E* > 3 GeV and in the polar angle range 0.35 < Q 1 < 
2.4 rad. The full acceptance involves all the event selec- 
tion criteria. 

As calculated using the AfkQed generator and the full 
simulation, the acceptance needs to be corrected for the 
effects resulting from the imperfect description of NLO 
ISR. For this correction Phokhara and AfkQed are com- 
pared at the generator level. Since the effect of the NLO 
differences is to give different longitudinal boosts to the 
events, one expects deviations in the geometrical and mo- 
mentum acceptance. This justifies the use of the genera- 
tors at 4-vector level. To improve on this, track and pho- 
ton parameters are smeared using resolution functions 
obtained from data. The acceptance is defined at this 
level by the polar angle ranges for the ISR photon and 
the two muons, and the p > 1 GeV/c requirement on the 
muons. 

We test the sensitivity of the results to using only 
fast simulation. Smearing generates a relative shift of 
1.0 x 10~ 3 for the acceptance correction. So any inad- 
equacy of the resolution functions is expected to be at 
a lower level. Some effects are not included in the fast 
simulation while they enter the full MC. The main com- 
ponents of the overall efficiency for the full simulation 
are shown in Fig. [27] While the loss of acceptance is es- 
timated to be 92% near the p mass, the fast simulation 
accounts for 78% only. The major contribution to the dif- 
ference between full and fast simulation comes from the 
DIRC crack removal and the IFR active area (cracks and 
bad areas). These azimuthal effects are essentially in- 
sensitive to the longitudinal boosts from additional ISR 
photons. Fast simulation is consequently adequate to 
describe the event acceptance changes generated by the 
additional ISR photon kinematics. 
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FIG. 27: Breakdown of the full simulation acceptance with 
respect to the generated events in the ISR photon angular 
range in the e + e~ center-of-mass 20° — 160° for p/u(7)7isr 
events. The numbers refer to the sequential application of 
the selection requirements: (1) trigger + acceptance selection 
for reconstructed ISR photon and tracks; (2) preselection of 
ISR events + E* > 3 GeV; (3) p > 1 GeV/c for both tracks; 
(4) tracks in IFR active area; (5) tracks in DIRC active area; 
(6) '/n/n'-ID + x 2 selection + J/ip rejection + minor selections. 



4- Effect of collinear additional ISR in AfkQed 



The angular distribution of hard additional ISR pho- 
tons can produce a significant transverse momentum that 
affects the event acceptance and preselection efficiency. 
The change of acceptance for collinear and non-collinear 
additional ISR is investigated with Phokhara as it pro- 
vides an MC sample with additional ISR following the 
QED angular distribution. A significant decrease of the 
acceptance is observed as a function of the polar angle 
of the additional hard (> 0.2 GeV) ISR photon. The dif- 
ference between Phokhara and AfkQed is aggravated by 
the nix-yi SIi > 8 GeV/c 2 requirement used at generation in 
AfkQed, which suppresses hard additional ISR. Both of 
these effects are kinematic in nature, and are well studied 
at the 4-vector level. 

The observed differences between data and AfkQed 
for the angular and energy distributions for NLO ISR 
(Sect. IV A2[) are overcome in Phokhara, which provides 
a much better description of the data. The geometri- 
cal acceptances computed with the smeared 4-vectors in 
Phokhara and AfkQed differ by about 2% in most of the 
mass range, Phokhara leading understandably to a lower 
acceptance (Fig. [25]). The global efficiency e MM ( 7 ) 7/SH ob- 
tained with AfkQed is corrected by this factor when com- 
puting the /X/u(7)7isr cross section. 
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FIG. 28: The ratio of the /J./-i(-y)~/iSR acceptances determined 
in Phokhara and AfkQed at the generator level with fast sim- 
ulation. 



B. ISR photon efficiency correction 

A coarse 7isr detection efficiency map in the (E*, 7 ) 
plane is derived from /i + /i - 7isR events selected on the 
basis of the muons only in data and full MC simulation. 
The data/MC correction for the ISR photon efficiency is 
obtained as a function of m MM by sampling the efficiency 
maps using the simulated sample. The efficiency is found 
to be lower in data by (1.5 ± 0.1)% below 2 GeV/c 2 , with 
a difference slightly smaller above. A systematic uncer- 
tainty of 3 x 10~ 3 is assigned to cover the effects origi- 
nating from the limited map granularity. 



C. Distributions of kinematic variables 

The comparison of distributions of relevant kinematic 
variables (polar angle 6 1 of the ISR photon, angular and 
momentum distributions of the muons) observed in data 
and simulation is an important cross-check of the analy- 
sis, as the true distributions are predicted by QED. Not 
all detailed corrections that are applied to the full simu- 
lation as a function of the \i\x mass are available for these 
variables, and we only consider corrections from PID for 
this test. Knowing some deficiencies of AfkQed for ad- 
ditional ISR, the comparison is made for events without 
excessive extra radiation, requiring the x 2 of the 1C fit 
that uses only the two tracks to be less than 15. 

Figures [29] and [30] show the distributions of B 1 and 
Pft in three mass intervals. In each case, MC is normal- 
ized to data as we are interested in testing the shapes. 
The agreement between data and the simulation is good, 



except for the ISR photon distribution at small an- 
gles where the data lies below the simulation. This ef- 
fect, which cancels in the nn/fifi ratio, is imputed to a 
data/MC difference for the photon efficiency at small an- 
gles. 

The angular distribution (0* ) of the muons in the \i\x 
cm. with respect to the ISR photon direction in this 
frame is of particular interest since it is predicted by QED 
to behave as 



dN 
dcosO* 



■ 1 



- cos 2 



£ + (l-/3 2 )sm 2 0* , 



(15) 



for pure ISR production, with the muon velocity /3 = 

1 — 4m 2 / s' . So we expect the distribution to be flat 

at threshold and to follow a 1 + cos 2 0* distribution at 
intermediate mass. At higher masses a larger fraction 
of the 'ISR'-selected photon comes in fact from FSR, in- 
creasingly modifying the cos 0* distribution. 

The distributions of | cos #* | for different mass intervals 
agree well with expectation as seen in Fig. [31] Although 
they are strongly biased by the p > 1 GeV/c requirement, 
which depletes the region near one, distributions in the 
threshold region indeed show the behavior expected from 
Eq. (fT"5|) , in agreement with the MC. 

Thus the distributions of the main kinematic variables 
of the selected ^^(^)^isr sample are in good agreement 
with expectations from QED as implemented in the MC 
generator. 



D. Systematic uncertainties on the absolute 
w(7)7ish. cross section 

The statistical errors of the measured efficiencies are 
included with the main statistical uncertainty on the /x/i 
mass spectrum. However, in some cases, remaining sys- 
tematic uncertainties are attached to the efficiency mea- 
surement process. Estimated systematic uncertainties 
on the measured cross section are summarized in Ta- 
ble [IT] for the mass range from threshold to 2GeV/c 2 . 
Above 2 GeV/c 2 the uncertainties are smaller, essentially 
because of the more straightforward determination of the 
muon-ID efficiencies. In some cases no systematic error 
is quoted when all uncertainties proceed from measure- 
ments and are already included in the point-to-point sta- 
tistical errors. 

We find that Phokhara at the fast simulation level is 
adequate to correct the AfkQed generator and that the 
systematic uncertainty on the acceptance resulting from 
the fast simulation is 10~ 3 . The effect of the momentum 
calibration uncertainty is only at the 10 -4 level. 

The absolute /i/i (7)tisr cross section makes use of the 
effective luminosity function defined by Eq. ^ that in- 
cludes the BABAR luminosity L ee . The latter is obtained 
for all the analyzed data using measurements of ee — > ee, 
/i/i and 77, and amounts to 230.8 fb^ 1 . The correspond- 
ing systematic uncertainty is 0.94%. 
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FIG. 29: The comparison between the distributions of data 
(points with errors) and simulation corrected for data/MC 
differences in PID (blue histogram), for # 7 in radians in the 
m MM intervals 0.5 - 1 GeV/c 2 (top), 1.5 - 2GeV/c 2 (middle), 
2.5-3GeV/c 2 (bottom). 



FIG. 30: The comparison between the distributions of data 
(points with errors) and simulation corrected for data/MC 
differences in PID (blue histogram), for ±p M ± in GeV/c in the 
m MM intervals 0.5 - 1 GeV/c 2 (top), 1.5 - 2 GeV/c 2 (middle), 
2.5 - 3GeV/c 2 (bottom). 
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TABLE II: Systematic uncertainties (in 10 ) on the ab- 
solute /^(7fsr) cross section from the determination of the 
various efficiencies in the fifj, mass range up to 2 GeV/ c 2 . The 
statistical part of the efficiency measurements is included in 
the total statistical error in each mass bin. For those con- 
tributions marked '-' all the relevant uncertainties come from 
measurements and are already counted in the statistical er- 
rors. 



Sources 


systematic errors (10 3 ) 


triggers and background filter 


0.3 


tracking 


1.3 


muon ID 


3.3 


tttv and KK backgrounds 




multihadronic background 




\ 2 cut efficiency 




angle and momentum acceptance 


1.0 


ISR photon efficiency 


3.4 


e + e~ luminosity 


9.4 


NNLO corrections to ctqed 


2.0 


sum 


10.9 



icose'i 



The uncertainty assigned to the QED cross section 
comes from the neglect of the NNLO contribution. The 
latter is estimated from the NLO fraction as given by 
PHOKHARA, equal to (4.33 ± 0.11)% within the selec- 
tion used in this analysis. Assuming a geometric growth 
of the coefficients of the expansion in a, the NNLO frac- 
tion is estimated to be 2 x 10~ 3 , which is taken as the 
systematic uncertainty. 

Summarizing, the overall systematic uncertainty on 
the absolute ^^(-j)jisr cross section is 1.1%, dominated 
by the BABAR luminosity error. 




lcos6*l 



FIG. 31: The comparison between the distributions of data 
(points with errors) and simulation corrected for data/MC 
differences in PID (blue histogram), for Icos^'l in the 
intervals: (from top to bottom) 0.20 - 0.25 GeV/ c 2 , 0.25 - 
0.30 GeV/c 2 , and 0.5 - 1 GeV/c 2 . 



E. QED test with the fifi^yisn events 

The comparison of the hh(^)^isr cross section with 
QED is made through the ratio of the distribution in 
the data as a function of to the same distribution 
of the simulation. Specifically, the distribution of the 
data is background-subtracted, and the distribution of 
the AfkQed-based full simulation, normalized to the data 
luminosity, is corrected for all data/MC detector and re- 
construction effects and for the generator NLO limita- 
tions using the Phokhara/AfkQed comparison with fast 
simulation. Because of the latter adjustments, discussed 
in detail in Section IVIII A[ the corrected ratio of m PM 
spectra is equivalent to a direct comparison of data with 
the NLO QED cross section. 

The QED prediction for the m = distribution is 
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obtained in the following way: 



cLNqed 

dm 



T ^r NL 

J-'ee Ophokhara 



1 dN\ 
N dm J 



AfkQcd,M>8 



fullsim 



J_dN_\ 



{n 



Phokhara 



dm J f + • 

, v AfkQcd m M>8 X Cdata/MC , (16) 

(j_dN_) 

where for each case Nq is the generated number of events, 
dN I dm the mass spectrum of events satisfying all cri- 
teria. The ratio of spectra at the generator level with 
fast simulation are labeled 'fastsim', while 'fullsim' de- 
notes the spectrum of events with full detector simula- 
tion. The notation 'M > 8' recalls that AfkQed was 
run with a requirement limiting hard additional ISR, 
tox 7isr > 8 GeV/c 2 . Finally, the Cd a ta/MC factor incor- 
porates all data/MC corrections for detector efficiencies, 
such as trigger, tracking, muon ID, \ 2 selection, and ISR 
photon efficiency. 

Because the PID efficiency varied with time due to the 
degradation of the IFR detector, the ratio data/QED 
is determined separately splitting the running period 
in two. Both distributions are flat from threshold to 
3.5 GeV/c 2 and consistent with unity within errors with 
satisfactory x 2 values. The difference of the ratios for 
the two data sets is (6.0 ± 4.0 ± 3.5 ± 4.4) x 10" 3 , where 
the first error is statistical, the second from non-common 
systematics (uncorrelated parts of the //-ID systematic 
uncertainties), and the third from the BABAR luminos- 
ity. Since the two samples correspond to different per- 
formances of the IFR detector, this test provides a con- 
firmation that the muon-ID efficiency has been handled 
adequately. 

The fit to the full data set is shown in Fig. [33] and 
yields 



G 



data 

W(7)7ISR 
NLO QED 
MM(7)7isr 



1 + (4.0 ± 1.9 ± 5.5 ± 9.4) x 10~ 3 , (17) 



where the first error is statistical, the second from sys- 
tematics, and the third from the BABAR luminosity. The 
value found for the ratio is consistent with unity over the 
full mass range explored in this analysis. We conclude 
that the measurement of the e+e~ — > ^ (7)7isr cross 
section using the BABAR luminosity agrees with NLO 
QED in the /z/x mass range from threshold to 3.5 GeV/c 2 
within an overall accuracy of 1.1%. 



F. Determination of ISR luminosity 

In this section we express the results obtained for the 
A t A t (7)7iSR sample in terms of the effective ISR luminos- 
ity, following Eq. ([TJ . As discussed in Section IVIII Al 
the jU/x(7)7isr event acceptance, ew7)7 JSJi > appearing 
in Eq. (JlJ is obtained from the large simulated sample 
generated with AfkQed with corrections for detector and 
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FIG. 32: The ratio of the \x\i mass spectrum in data over the 
absolute prediction from QED using the BABAR luminosity. 
The NLO QED prediction is obtained from the data-corrected 
(for detector simulation) and Phokhara-corrected (for NLO 
effects) AfkQed mass spectrum. The band is drawn around 
the fit of the 0.2-3.5 GeV/c 2 region to a free constant, with a 
width given by ± the total expected systematic uncertainty. 



reconstruction effects. Corrections specific to the QED 
test, i.e., NLO and ISR photon efficiency corrections that 
cancel in the tttt/ /i/j, ratio, are not applied. 

Several factors need to be considered in addition: (i) 
the LO FSR correction, (ii) unfolding of the data from 
to v / s' to include the possible emission of an ad- 
ditional FSR photon, and (iii) the QED cross section 



M + A* (7fsr)) ( s ') at the Born level concern- 



ing ISR, but including FSR. Except for (ii), which has 
been discussed in Section [VII Al we address these points 
in turn before giving the final result on the ISR luminos- 
ity. 



1. Lowest-order FSR correction 

The most energetic detected photon is assumed to be 
emitted by the initial state. This is largely true at low 
mass, but there is an increasing probability at larger 
s' values that this photon originates from muon radia- 
tion. Thus the observed /x/z mass spectrum has to be 
corrected in order to keep only ISR production, since for 
all practical purposes at BABAR where t/s ~ 10.58 GeV 
and Vs^ < 5 GeV, LO FSR production is negligible for 
hadronic processes. 

Figure [33J shows the <5pg R correction to the cross sec- 
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FIG. 33: The FSR correction 5^ R obtained with AfkQed. 



tion, defined as 



°FSR 



1-4 



FSR 



.4 



add.ISR.add.FSRl 



1-4 



ISR 



A 



add.ISR.add.FSRl 



(18) 



as a function of ys', where .Afsr (-^isr) is the LO FSR 
(ISR) amplitude and -4 a dd.iSR,add.FSR is the NLO contri- 
bution. <5pg R is obtained with AfkQed at the generator 
level. It would be preferable to use Phokhara instead, 
as we know additional ISR generation is approximate in 
AfkQed, but by construction the FSR or ISR origin of 
photons is not available in Phokhara, hence s' is not ac- 
cessible on an event-by-event basis. However the differ- 
ence in <5pq D is expected to be at a negligible level, about 



J FSR 

and 2 x 



10 6 at 1 GeV and 3 GeV, respectively. 



2. Born QED cross section with additional FSR 

The cross section for e + e~ — > (7fsr); at the Born 

level for the initial state and without vacuum polariza- 
tion, can be calculated exactly in QED at NLO. It has 
the form: 



with 



= (1 + Cd.FSR) 

47ra 2 (0) 13(3 ~/3 2 ) 



o-pt(s') = 



3s' 



"add. FSR 



a(0) 



V( a> ) = Vh(s') + Vs(s') + Vv(s') , 



(19) 



(20) 

(21) 
(22) 



where /3 is the muon velocity and rjhsv are the 
0(a) contributions to the final state from hard and 



soft bremsstrahlung, and the one- loop/Born interference 
('virtual' contribution), respectively. 

The sum of r\ v and r\ s is infrared-finite, while the total 
sum is independent of the choice of the energy used to 
separate soft and hard photons (within reasonable lim- 
its). Expressions for all three components can be found 
in many papers, for example in Refs. [25l [26j. Byvirtue 
of the Kinoshita-Lee-Nauenberg (KLN) theorem 27j, the 
dominant logarithmic terms cancel in the sum of the 
(soft+virtual) and hard contributions. Although the two 
terms reach a level of a few percent, they have opposite 
signs and the sum S^ d FSR stays in the few x 10~ 3 range. 
This explains why a sizeable hard additional-FSR signal 
is seen in data, despite the fact that the total additional- 
FSR contribution is very small. 



3. The effective ISR luminosity for the 7T7t(7)7isr analysis 

For the 7T7r(7)7isR analysis, the luminosity dLfg R inte- 
grates all configurations up to two ISR photons, where 
at least one photon has E* > 3 GeV and is in the an- 

c.m. with 6>* lin = 



mm ' 

= 20° 



gular range 

180° - 6»* lax _ 

The full effective ISR luminosity dL^ K /dy/s' is de- 
rived from the measured /^//(7)tisr spectrum according 
to Eq. ([!}. The event acceptance is taken from AfkQed, 
with corrections for detector and reconstruction effects. 
Unfolding of the background-subtracted /i/i mass spec- 
trum is performed as explained in Section IVII A[ The 
result is shown in Fig. [3U The effective luminosity de- 
rived this way implicitly includes the VP factor since 
the /x/z(7)tisr data include vacuum polarization effects, 
while the bare cross section tT ° M ( 7 )(v / s') entering Eq. ([1]) 
does not. 

The measured effective luminosity is compared to the 
standard estimate of Eq. ([B]) using LO QED, given by 



dL 



LO 



dW LO fa(s') 



dV? Lec dV? W«) 

where the LO radiator function is [l^, [ll[ 
dW LO a(0) 



(23) 



7r.r 



(2 - 2a; + ar) In ■ 



with x = 1 — s'/s and c = cos0j^ in . We insert the VP 
term a(s')/a(0) in Eq. (|23[) for a convenient comparison 
with the effective luminosity. The computed VP factor 
includes both leptonic and hadronic contributions. The 
hadronic contribution is taken from the parametrization 
used in AfkQed and is found to agree well with an inde- 
pendent determination using the tools of Ref. 0. 

The LO+VP prediction is superimposed on the mea- 
sured effective luminosity in Fig. 1341 The measured lumi- 
nosity is found to be about 2% larger than the LO+VP 
QED result. This difference varies slowly with mass and 
includes systematic effects on the BABAR luminosity de- 
termination, the effect of the NLO contribution in data, 



39 



4000 



2000 





ii, 


ii 




• eiierave lumi. 
- lumi.(LO) + VP 
























12 3 

n/F ( GeV) 

FIG. 34: The effective ISR luminosity for the rnr analysis: the 
data points give Lisr in A-v/s' = 50 MeV bins. The conditions 
for the detected/identified ISR photon are Ej > 3 GeV and 
20° < 61 < 160° in the e + e~ cm. frame, while one additional 
ISR photon is allowed without any restriction. The superim- 
posed histogram is the lowest-order ISR prediction following 
Eq. (|23p . The J/ip mass region is removed for the luminosity 
determination. 



the difference between the ISR photon efficiency in data 
and MC and any residual effect in the detection efficiency. 
The latter contribution is small, in accordance with the 
successful QED test performed with the W.( 7 ) 7 isr cross 
section. 

The effective luminosity shown in Fig. [M]is measured 
in 50 MeV bins. This interval size is too wide near nar- 
row resonances (w and 4>) because of the rapid variation 
of hadronic vacuum polarization. Therefore, we com- 
pute the local variation inside each 50 MeV interval as 
the product of the lowest-order QED radiator function 
times the vacuum polarization factor. The result is nor- 
malized to the effective luminosity in the interval. In this 
way the detailed local features are described, while pre- 
serving the measured effective luminosity as a function 
of mass. 

The statistical error of the irir cross section is limited in 
the p resonance region by the number of events available 
to determine the ISR luminosity. Bin-to-bin statistical 
fluctuations are reduced by a suitable averaging of the 
ratio of measured to LO ISR luminosities. The ratio dis- 
tribution in 50 MeV bins is smoothed by averaging five 
consecutive bins in a sliding way. This value is chosen 
as a compromise between smoothing and the validity of 
the assumption of slow variation. This method does not 
improve in principle the ISR luminosity statistical error 
because the reduced local error is compensated by the 
correlation between neighbouring bins. A slight improve- 
ment in the dispersion integral is however observed due 



to the weighting in different mass regions. 

The statistical errors on the ISR effective luminos- 
ity from the measurement of efficiencies are included in 
the statistical covariance matrix, while the systematic 
errors from the different procedures are accounted for 
separately. These errors are 0.3 x 10~ 3 for the trig- 
ger, 1.3 x 10~ 3 for tracking, 2.9 x 10~ 3 for //-ID, and 
1.0 x 10~ 3 for acceptance, for a total systematic uncer- 
tainty of 3.4 x I0 -3 . The uncertainty from the correlated 
loss of /z-ID for both tracks is not included here, since it 
is anticorrelated with the pion rate. It is counted in the 
systematic errors on the 7T7t cross section. 



IX. MEASUREMENT OF THE e+e" ->■ 7r+7r-( 7 ) 
CROSS SECTION 

The e + e~ — > 7r + 7r~(7) bare cross section is measured 
from the y/s' distribution of produced tt + tt~ (7)7isr 
events divided by the effective ISR luminosity obtained 
from /x + /i~( 7 ) 7 igR on the same data. The v / s' distribu- 
tion is obtained from the observed m^^ mass spectrum, 
after background subtraction, corrections for data/MC 
efficiency differences, unfolding, and MC acceptance cor- 
rections, as described in detail in the preceding sections. 
Because the cross section spans several orders of magni- 
tude over the energy region considered, from threshold 
to 3 GeV, the analysis strategy depends on the mass re- 
gion. Event selection and background subtraction are 
optimized separately for the p resonance central region 
or for the resonance tails, with corresponding efficiency 
corrections and unfolding matrix. To facilitate compar- 
ison with other experiments, the results are also shown 
in terms of the pion form factor fitted with a vector- 
dominance model (VDM). 

A. The central p region (0.5 < m nn < lGeV/c 2 ) 
1 . Strategy 

The mass region between 0.5 and 1 GeV/c 2 , dominated 
by the p resonance, provides the dominant contribution 
to the vacuum-polarization dispersion integrals and their 
errors. The need for small systematic uncertainties, con- 
gruous with the small statistical errors, together with the 
low background level in that region, lead to an event se- 
lection with the largest efficiency. Therefore the loose \ 2 
criterion, the same as for the PP^Jisr analysis, and 
standard 7T-ID for both tracks are used. 



2. Summary of backgrounds 

The backgrounds are obtained as described in Sec- 
tion IVII The dominant contribution is from multi- 
hadronic processes, mostly ISR (7r + 7r~7r°7, 7r + 7r _ 27r°7) 
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and qq, with a fraction amounting to 8.4 x 10 3 at the p 
peak. The ppj contribution is much smaller (< 10 -3 ). 

TABLE III: Estimated background fractions (in %) in the 
'tttt' sample for 7^=0.525, 0.775, 0.975 GeV/c 2 . The quoted 
errors include both statistical and systematic uncertainties. 



Process 0.525 GeV/c 2 0.775 GeV/c 2 0.975 GeV/c 2 





3.48 ± 0.36 


0.37 ±0.23 


2.71 ±0.31 


KK 


0.08 ±0.01 


0.01 ±0.01 


0.08 ±0.01 


727T7T 


8.04 ±0.41 


0.39 ±0.05 


0.88 ±0.19 


qq 


1.11 ±0.17 


0.26 ±0.03 


1.81 ±0.19 


727T27T 


1.29 ±0.16 


0.06 ±0.01 


0.46 ±0.09 


747T 


0.20 ± 0.04 


0.09 ±0.01 


0.24 ±0.06 


7PP 


0.22 ±0.02 


0.04 ±0.01 


0.52 ±0.06 


yrfL-K 


0.02 ±0.01 


0.03 ±0.01 


0.09 ±0.01 


~jK s K L 


0.18 ±0.03 


0.01 ±0.01 


0.10 ±0.02 


747T27T 


< 0.01 


< 0.01 


< 0.01 


TT 


0.17 ±0.03 


0.04 ±0.01 


0.31 ±0.05 


7ee 


0.63 ±0.63 


0.03 ±0.03 


0.27 ±0.27 


total 


15.38 ±0.87 


1.31 ±0.24 


7.37 ±0.51 
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FIG. 35: The total background fraction for the 7T7r(7)7isR 
sample in the central p region (loose \ 2 selection). 



The fractions of all the considered backgrounds are 
given in Table Hill at three mass values. For convenience, 
we also show the level of the pipi and KK background con- 
tributions in the '7r7r'-idcntificd sample, although they 
are implicitly subtracted when solving Eqs. (j8J) for the 
produced spectrum N n7T . 

The total background fraction as a function of ra^n is 
shown in Fig. |35] It is 1.3% at the p peak, but reaches 
-15% at 0.5GeV/c 2 and -7% at lGeV/c 2 . These sharp 
increases justify the limits chosen to define the 'central 
region'. At the worst place, near 0.5GeV/c 2 , the total 
uncertainty from the estimated non-fip/KK background 
fraction is 0.8%, which is still tolerable. At the peak the 
uncertainty is less than 0.1%. 
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3. Background-subtracted m„ mass distribution 

The background-subtracted m n7r spectrum obtained in 
the p region before unfolding, with loose % 2 criterion and 
l 7r'-ID for both pions, is shown in Fig. (35] Only the sta- 
tistical errors in the 2MeV/c 2 mass intervals are given, 
amounting to 1.4% on peak and 4.4% near the bound- 
aries. Apart from the p resonance shape, a clear p — lj 
interference pattern is observed. 



B. The p tail regions {m^Tt < 0.5, m W7r > lGeV/c 2 ) 
1 . Strategy 

The pion cross section decreases very rapidly away 
from the p resonance, while the backgrounds from 
wilhiSR, KK(j)ji SR , ppjisr, and multihadron 



0.5 



0.6 



0.7 



0.8 



0.9 1 
m^GeV/c 2 ) 



FIG. 36: The 
region, in 2 MeV/c 2 I 



spectrum of 7T7r(7)7isR events in the p 



events show a smooth variation with the 7T7r mass. To 
keep the background levels under control in the p tail 
regions, the selection of ttttI^jisr events is tightened 
with respect to the criteria used in the p region. Two 
handles are simultaneously used: (i) the tight \ 2 crite- 
rion ln(x 2 gR ±1) < 3 is chosen to reduce multihadronic 
backgrounds, and (ii) the pion-ID is strengthened to im- 
prove muon (and also electron) rejection. In addition, 
the Vxy requirement described in Section IVI B 41 is ap- 
plied to remove backgrounds from photon conversions 
and bremsstrahlung in the beam pipe. 
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2. Summary of backgrounds 

For two-body ISR and ee backgrounds, the tighter x 2 
criterion is not useful, so the harder tt identification, 
'iTh (Sect. IIVD 2[) is required for at least one of the 
two V-identified tracks, giving a further rejection of p 
and e. The downside is that, because the V^' identifica- 
tion breaks the completeness of PID classes, the pp, and 
KK backgrounds cannot be subtracted anymore from 
the l 7r7r^' sample by solving the Eqs. (0 system. 

The pp{^)^fisR background is estimated according to 
Section IVlB 11 With the 'irirh' selection, it is reduced by 
a factor ~ 7 with respect to the background remaining 
after the standard pion-ID. However the corresponding 
factor for the background fraction in the final pion sample 
is smaller because the 'mr^ efficiency also reduces the 
signal (Sect. IIVD 711 . 

Neither the tight \ 2 criterion nor the 'with' ID brings 
significant reduction of the KK(j)jisr background com- 
pared to the selection used in the p peak region. Its small 
contribution is estimated from data by using the proce- 
dure described in Section IVl B 21 and subtracted. 

The backgrounds from qq and multihadronic ISR 
events are estimated as discussed in Sections IVI B 61 and 
IVI B 71 respectively. They are much reduced compared 
to the central p region because of the tight \ 2 condition. 

The different fractions of background in the region of 
the p tails, with '7r7iV ID and hi(x 2 SR + 1) < 3, are given 
in Fig. [37l Fractions at specified masses are listed in Ta- 
ble QV] The total background contribution is obtained by 
summing all the individual contributions obtained above. 

3. Background-subtracted m-irir mass distribution 

The background-subtracted m n7r distribution of 
7T7t(7)7isr events before unfolding, using 'tttt;,' identifi- 
cation and ln(xFSR + -0 < 3 is plotted from threshold to 
3 GeV/c 2 in 50 MeV/c 2 mass intervals in Fig. [33 A dy- 
namic range of 10 3 -10 4 is observed between the p peak 
and either the first bin above threshold or at 3GeV/c 2 . 
The dip structure at 1.6 GeV/c 2 seen by the DM2 exper- 
iment [281 ] is confirmed with high statistics and a new 
structure shows up near 2.2GeV/c 2 . 

C. Angular distribution in the nn center-of-mass 

The distributions of kinematic variables such as the 
ISR photon polar angle and the pion momenta and angles 
depend on the hadronic structure we seek to measure. 
Thus the detailed comparisons between data and MC 
distributions expected from theory, which are performed 
in the ppirfjisR channel with QED (Sect. IVIII Cp . are 
meaningless in the pion channel. However, one distri- 
bution, namely the pion angular distribution in the tttt 
center-of-mass, with respect to the ISR photon direc- 
tion in that frame, is model-independent. The cos 6* 



distribution is consistent with sin 2 #* as expected in the 
e + e _ — > 7r + 7r~ process, but it is strongly distorted at 
I cos 9* I values near one by the p > 1 GcV/ c requirement 
on the tracks. 

The I cos #* I distributions for background-subtracted 
data and MC are compared in Fig. [39] for the 0.5- 
1 GeV/c 2 mass range: they agree with each other within 
the statistical errors, as expected for a pure pion sample. 



D. Acceptance and corrections 

The overall efficiency e^Mtisr entering Eq. (0 for the 
pion channel is calculated using the AfkQed generator 
and full simulation in the same way as for the muon chan- 
nel. It is corrected for differences in efficiencies between 
data and MC, which are introduced as mass-dependent 
corrections applied to the event spectrum (Eq. (0). 

As discussed in Section IVIII Al NLO approximations 
are made in simulation, which affect the acceptance. 
While the FSR prescription by PHOTOS is found to agree 
reasonably well with data, this is not the case for addi- 
tional ISR as simulated in AfkQed. The problems have 
been studied in detail for muons, since they affect the ab- 
solute measurement of the A*m(7)7isr cross section and 
the comparison with QED. However here we deal with ac- 
ceptance corrections that apply to the pion cross section 
measured from the ratio of the 7T7t(7)7isr spectrum to 
the effective luminosity. As the additional-ISR issues are 
common to the 7T7t(7)7isr and nnirfyisn channels, they 
cancel in the i? oxp (v / s') ratio, except for second-order ef- 
fects addressed below. Thus the tttt measurement does 
not rely on the accurate description of NLO effects by the 
MC generator, a fact that is a strength of this analysis 
method. 

Acceptance, including preselection efficiency, is mainly 
affected by kinematics, i.e., the angular and energy dis- 
tributions of the hard additional-ISR photon. We study 
these effects for the pion channel as we have done for the 
muon channel, using AfkQed and Phokhara at the gener- 
ator level and fast simulation with parametrized efficien- 
cies and resolutions. The resulting data/MC correction 
on acceptance for pions is consistent with the correction 
obtained for muons. Second-order corrections induced by 
pion secondary interactions are investigated, using the 
full simulation. The total acceptance correction amounts 
to a few x 10~ 3 in the p mass region. It is shown as a 
function of the tttt mass in Fig. 001 As for the other cor- 
rections, the statistical uncertainties are included in the 
final cross section errors. Since part of the correction is 
derived using the fast simulation at the generator level, 
a significant fraction of the correction (~25%) is taken 
as a systematic uncertainty, 10 -3 in the 0.6-0.9 GeV/c 2 
region and larger outside. 
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FIG. 37: The fractions of different backgrounds in the physical sample with the tight 2D-x criterion and strengthened l 7T7iV-ID 
(as used in the p tails region) as a function of the -kit mass, (top loft): Multihadrons, including tt. (top right): /i/i(7)7isR 
(data + measured mis-ID), (bottom left): KK{^f)^fisn (data + measured mis-ID), (bottom right): pp7isr (MC). 



TABLE IV: Estimated background fractions (in %) in the l 7r7r h ' sample for 771,^=0.325, 0.475, 0.975, 1.375, 1.975, and 
2.975 GeV/c 2 . The entries marked as ' — ' correspond to a negligible fraction. Processes with fractions less than 0.05% in all 
intervals are not listed. The quoted errors include both statistical and systematic uncertainties. 



Process 0.325 GcV/c 2 0.475 GcV/c 2 0.975 GeV/c 2 1.375 GeV/c 2 1.975 GeV/c 2 2.975 GeV/c 2 



fifj, 


7.7 ±2.5 


1.4 ±0.4 


0.5 ±0.1 


3.6 ±0.6 


10.5 ± 1.9 


56.2 ± 15.8 


KK 


8.2 ±0.7 


0.3 ±0.1 


0.1 ±0.1 


0.4 ±0.1 


0.2 ±0.1 


0.2 ±0.2 


727T7T 


0.4 ±0.2 


0.3 ±0.1 


0.1 ±0.1 








qq 


0.1 ±0.1 


0.1 ±0.1 


0.3 ±0.2 


5.0 ±2.2 


0.8 ±0.7 


3.4 ±4.5 


1PP 


0.7 ±0.1 


0.3 ±0.1 


0.4 ±0.1 


4.3 ±0.5 


9.2 ± 1.4 


24.5 ±6.8 


total 


17.1 ±2.6 


2.3 ±0.4 


1.4 ± 0.2 


13.2 ±2.3 


20.7 ±2.5 


84.3 ± 17.8 
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FIG. 38: The m^vr spectrum of 7T7r(7)7isR events selected with 
'iTTTh' identification and the tight ln(xisR + 1) < 3 criterion, 
from threshold to 3 GeV/c 2 in 50 MeV/c 2 mass intervals. 




Icos0*l 

FIG. 39: The angular pion distribution in the nn system with 
respect to the ISR photon direction as function of j cos 8% \ 
for background-subtracted 7T7t(7)7isr data (points) in the p 
central region (0.5 < < lGeV/c 2 ). The blue histo gram is 
the shape obtained in the simulation, normalized to the data. 



E. Summary of the treatment of statistical 
uncertainties 

The statistical covariance matrix of the cross section 
includes the bin-to-bin correlations affecting the 7T7t spec- 
trum and the luminosity. 

The statistical covariance matrix of the 7T7r spectrum is 
not diagonal, due firstly to correlations introduced by the 



FIG. 40: The full correction to the £nn(-y)j ISR /E f j,v(i)yisR ac ~ 
ceptance ratio for non-canceling effects (see text). 

transfers of events in the unfolding process. In addition, 
the data/MC efficiency corrections and subtracted back- 
ground spectra are initially computed in 50 MeV/c 2 bins, 
but applied to spectra with 2 MeV/c 2 bins (in the central 
p region) or 10 MeV/c 2 bins (in the p tail regions) using 
splines. The resulting covariance matrix is obtained from 
a large series of toy experiments. 

The ratio of the measured luminosity to the LO lu- 
minosity including vacuum polarisation is initially com- 
puted in (almost) uncorrelated bins of 50 MeV. The pro- 
cedure of sliding bins, used for smoothing this distri- 
bution, introduces correlations between the final values 
(Sect. IVllIF 3[) . The luminosity errors for the final cross 
section (2 or 10 MeV) bins are 100% correlated within 
a 50 MeV bin, whereas additional correlations occur be- 
tween the 50 MeV bins because of the bin-sliding pro- 
cedure. Finally, the correlation effect from unfolding the 
pp spectrum is rather weak, but it is however propagated 
to the final correlation matrix. 



F. Systematic errors 

Systematic uncertainties affecting the 7T7t sample in dif- 
ferent mass regions are now summarized. The statistical 
errors of the measured efficiencies are included with the 
main statistical uncertainty on the irir mass spectrum. 
However, in some cases, remaining systematic uncertain- 
ties are attached to the efficiency measurement process 
and quoted as such. Details have been given for each ef- 
ficiency study in Sections IIVI and l"V Bl The results for all 
systematic uncertainties are listed in Table |V] 

The overall relative systematic uncertainty on the 
7T7t(7fsr) cross section is 5.0 x 10~ 3 in the 0.6-0.9 GeV 
range, but significantly larger below and above the cen- 
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FIG. 41: The ratio of the corrected and unfolded mass spec- 
tra (data points) for loose over tight 2D-x 2 selection in the 
central p region fitted in 100 MeV/c 2 bins, compared to the 
band of independently estimated uncertainties (solid lines). 
The MC mass-matrix resolution correction is shown as the 
dashed histogram. 

tral region. For comparison, the statistical error of the 
measured efficiency corrections amounts to4.7xl0~ 3 at 
the p peak, while the statistical error of the raw spectrum 
is 1.35% at that mass. 

A full treatment of the systematic uncertainties is im- 
plemented, using a covariance matrix. To achieve this we 
consider the individual systematic errors (for each source, 
as given m TableE} to be 100% correlated in all the mass 
bins. Then the total systematic covariance matrix is built 
as the sum of the covariance matrices corresponding to 
each individual systematic source. 

G. Consistency check with tight and loose \ 2 
selection 

The loose \ 2 criterion is used in the p central region, 
while the tight one is used in the tails where backgrounds 
are larger. However it is possible to compare the results 
obtained with the two selections in the central region. 
This provides a test of the % 2 selection efficiency and of 
the multihadronic background. The test is also sensitive 
to unfolding, as mass resolutions are different in different 
2D-x 2 regions. For this test, events are selected with the 
l p central' conditions, and with either the tight or loose 
X 2 criterion. 

The result of the test is expressed as the ratio of the 
efficiency-corrected and unfolded spectra for the loose 
over the tight \ 2 selections. The fitted value of this ratio 
over the full central range (0.5-1.0 GeV/c 2 ) is found to 



be consistent with unity within errors, 0.9983 ± 0.0049 
with a x 2 /DF of 53.6/49 for 10MeV/c 2 bins. Fits in 
100MeV/c 2 intervals, given in Fig. [JTJ do not show any 
significant trend for a resolution mismatch between data 
and corrected MC. Deviations from unity are at a much 
smaller level than the resolution correction applied to 
the MC in the intermediate region (Sect. IVII A[ shown 
by the dashed histogram). They are also within the 
range of estimated uncertainties between the two \ 2 con ~ 
ditions (background and \ 2 selection efficiencies). Wc 
thus conclude that the procedure used for correcting the 
MC mass-transfer matrix is consistent within the quoted 
systematic uncertainties. 

H. Cross section results 

1. Estimates of LO FSR in the 7T7t(7)7 process 

In Eq. © it is assumed that the contribution of LO 
FSR to the irir(-f)j cross section is negligible. This ap- 
proximation is supported by calculations using specific 
models. 

The simplest model uses an extrapolation of the pion 
form factor to the large value of s ~ 112 GeV 2 , and as- 
sumes point-like pions to compute LO FSR photon emis- 
sion, as for additional FSR. This procedure, questionable 
for large energies E* > 3 GeV, yields a very small rela- 
tive contribution from |^4fsr| 2 3 ^fsr ~ 10~ 7 - A more 
realistic model 29] considers radiation from quarks and 
recombination into a pion pair, with the parametriza- 
tion of the produced even-spin states based on Ref. [30| . 
In this case the contribution is at a few x 10~ 4 level for 
masses below 1 GeV. However the FSR rate could be 
enhanced on specific resonances that are not exp licitly 
taken into account in the model. Estimates [31( using 
a 7*7/2(1270) transition form factor evaluated in the 
asymptotic regime by perturbative QCD indicate a FSR 
contribution of about 0.9% on the /2(1270) resonance. 
Contributions from /o(980) and /o(1370) are expected 
to be lower. 

Finally, a direct test with BaBAR data has been per- 
formed with a measurement of charge asymmetry, which 
is proportional to the interference between LO ISR and 
FSR amplitudes. This work in progress, which will be 
published separately, yields results that do not exceed 
the estimates above. 

The estimated LO FSR contributions are at levels 
smaller than the quoted systematic uncertainties on the 
7^(7)7 cross section, much smaller actually for the p re- 
gion. No subtraction has been applied to the measured 
cross section. 



2. Results on the bare cross section with FSR 

The results for the e + e~ — > 7t + 7t~(7fsr) bare 
cross section including FSR, er° , N , are given in 
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TABLE V: Systematic uncertainties (in 10~ 3 ) on the cross section for e + e~ — > 7r7r(7FSR.) from the determination of the various 
efficiencies in different ttty mass ranges (inGeV/c 2 ). The statistical part of the efficiency measurements is included in the total 
statistical error in each mass bin. The last line gives the total systematic uncertainty on the ivn cross section, including the 
systematic error on the ISR luminosity from muons. 



Sources 


0.3-0.4 


0.4-0.5 


0.5-0.6 


0.6-0.9 


0.9-1.2 


1.2-1.4 


1.4-2.0 


2.0-3.0 


trigger / filter 


5.3 


2.7 


1.9 


1.0 


0.7 


0.6 


0.4 


0.4 


tracking 


3.8 


2.1 


2.1 


1.1 


1.7 


3.1 


3.1 


3.1 


TT-ID 


10.1 


2.5 


6.2 


2.4 


4.2 


10.1 


10.1 


10.1 


background 


3.5 


4.3 


5.2 


1.0 


3.0 


7.0 


12.0 


50.0 


acceptance 


1.6 


1.6 


1.0 


1.0 


1.6 


1.6 


1.6 


1.6 


kinematic fit (x 2 ) 


0.9 


0.9 


0.3 


0.3 


0.9 


0.9 


0.9 


0.9 


correl pp ID loss 


3.0 


2.0 


3.0 


1.3 


2.0 


3.0 


10.0 


10.0 


Tin / fifi non-cancel. 


2.7 


1.4 


1.6 


1.1 


1.3 


2.7 


5.1 


5.1 


unfolding 


1.0 


2.7 


2.7 


1.0 


1.3 


1.0 


1.0 


1.0 


ISR luminosity 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


sum (cross section) 


13.8 


8.1 


10.2 


5.0 


6.5 


13.9 


19.8 


52.4 



Figs. 021 |43l and [44] as a function of \fp . The cross 
section is dominated by the wide p resonance, with struc- 
tures at larger masses. The dip region near 1.6 GeV, 
usually interpreted as resulting from interference be- 
tween the pf and p" amplitudes, is mapped with a much 
increased precision compared to previous experiments. 
There is also an indication for a structure in the 2.2- 
2.25 GeV region, which could be due to a still higher-mass 
p 1 " vector meson. 

Files containing the cross section data and their covari- 
ance matrices are provided in the EPAPS repository [32j . 



with 



/ / \ a ™h)( s ') ( a(s') \ 2 



and (3^ the pion velocity. The FSR correction [2^, [26[ , 
^add fsr = a (0)/ 7r v{ s ')i decreases slowly with s' and 
amounts to 8.0 x 10~ 3 at the p mass. 



I. Pion form factor fits 

The square of the pion form factor is defined as usual 
by the ratio of the dressed cross section without FSR, 
divided by the lowest-order cross section for point-like 
spin charged particles. Thus, 



I^IV) = 



3s' 



7ra 2 (0)£3 



(24) 



A vector-dominance model (VDM) is used to fit the 
BABAR pion form factor. It is a way to interpret the 
observed structures beyond the p resonance in terms 
of higher-mass isovector vector mesons. The fit also 
provides a convenient means to interpolate through the 
BABAR data points in order to facilitate the comparison 
to other experiments. 

The VDM parameterization, including p — u interfer- 
ence, is given by: 



1 + C p > + C p n 



+ c 



which satisfies F ff (0) = 1. The amplitudes of the Breit- 
Wigner (BW) functions are complex: c w = 



J<t>u 



-p' 



and c 



pin 



The BW of the uj is taken as: 



The wide p, p', p" and p'" resonances are described by 
the Gounaris-Sakurai(GS) model [33|, which takes into 
account the variation of their width with energy: 



BW* s (s,m,T) 



m 2 — s — imT 



(27) BW^(s,m,T 



m 2 (l + d(m)T/m) 



s + /(s, to, r) — imT{s, m, V) 



,(28) 
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FIG. 42: The measured cross section for e + e~ — > -k + -k~(^) over the full mass range. Systematic and 
statistical uncertainties are shown, but based only on the diagonal elements of the covariance matrix 
(see text). 
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FIG. 43: The measured cross section for e + e~ — > tt + -k~(j) in the lower mass range. Systematic and 
statistical uncertainties are shown, but based only on the diagonal elements of the covariance matrix 
(see text). 



with 



TV T-\ T S { ^( S ) 



(29) 



where /^(s) = yl — Arn^js. In principle this energy de- 
pendence is justified only below 1 GeV, as 4-pion final 
states dominate at larger energies, but it is used for sim- 



plicity. Detailed studies of the high mass states cannot be 
performed only on the basis of pion form factor fits, and 
require complex coupled-channel analyses. Such studies 
are beyond the scope of this paper, but the present 2tt 
data constitute a very useful ingredient for them. 



The auxiliary functions used in the GS model are: 
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FIG. 44: The measured cross section for e + e~ — > -K + TT~(-y) in the central p region. Systematic and 
statistical uncertainties are shown, but based only on the diagonal elements of the covariance matrix 
(see text). 



ji \ 3 ml 
dym) = — - 



7r A; 2 (to 2 ) 



In 



to + 2k(m 2 ) 



2nk{m 2 ) irk 3 (m 2 ) 



(30) 



f(s,m,T) = 



Tm 2 
k 3 (m 2 



[k 2 (s)(h(s) - h(m 2 )) + (to 2 - s)k 2 (m 2 )h'(m 2 )] , 



(31) 



r 



where 



= 2_m i J V 7s + 2k( S ) 

7T y'S 



2TO7T 



(32) 
(33) 



and h'(s) is the derivative of h(s). 

The form factor data is fitted in the full energy range, 
from 0.3 to 3.0 GeV, involving 18 free parameters: the 
mass and width of the p, and for each other resonance 
(oj, p', p", p'") the amplitude (modulus and phase) with 
respect to the p, and mass and width. According to a 
well-known effect [34l ] , the x 2 minimization returns fitted 
values that are systematically shifted with respect to the 
data points when the full covariance matrix is used in 
the fit, while the fit using diagonal errors is verified to 
be bias-free. This feature is due to correlations, which 
here arise from both statistical and systematic origins, 
but mostly from the ISR-luminosity 50 McV sliding bins 
(Sect. lViTIF 3[) and systematic errors. To circumvent the 
problem, we fit the data with only diagonal errors to ob- 
tain the central values of the fitted parameters. The error 



on each parameter is taken as the largest error obtained 
from the fit either with the full covariance matrix or with 
only diagonal errors. The biases on the mass scale cali- 
bration and the resolution obtained in Sections I VII Bl and 
IVH CI are included in the fit results on the p and ui res- 
onance parameters in Table IVH with the corresponding 
systematic uncertainties indicated. 

As shown in Fig. [JSJ the VDM fit provides an adequate 
description of the BABAR data over the full 0.3-3 GeV 
range (x 2 /DF = 351/319). The goodness of the fit shows 
that the GS parametrization of the dominant p resonance 
describes the data in a reasonable manner, as well as the 
contributions from the higher p', p" and p'" resonances. 
In particular the strong interference dip near 1.6 GeV is 
well reproduced. Beyond 2 GeV, the p'" is required in 
order to reproduce the structure seen in the data. The 
quality of the fit is shown in more detail in Fig. 05] in the 
low-mass range and in the p peak region with the p — uj 
interference. 

The relative ratio (I^Trldata/l-^rlvDM — 1) is shown in 
Fig. [IT] over the full energy range. Some deviation is 
observed in the low-mass region where the fit underes- 
timates the data. Some oscillation is also observed be- 
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FIG. 46: The pion form factor-squared measured by BABAR as a function of VV and the 
VDM fit from 0.3 to 3GeV described in the text, (top): Low-mass region (0.3-0.5 GeV). 
(bottom): p peak region with p — to interference (0.70-0.82 GeV). 
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FIG. 47: The relative difference between the pion form factor-squared from BABAR data 
and the 18-parameter phenomenological fit in three mass regions. Systematic and statis- 
tical uncertainties are included for data (diagonal errors). The width of the band shows 
the propagation of statistical errors from the fit and the quoted systematic uncertainties, 
added quadratically. 
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TABLE VI: Parameters obtained for the VDM fit (described 
in the text) to the BABAR pion form factor data. The errors 
include both statistical and systematic uncertainties. The 
errors shown in parentheses for the p and u> parameters stem 
from the mass calibration and resolution uncertainties (see 
text). 



Parameter 


Value ± Error 


m p (MeV/c 2 ) 


775.02 ± 0.31 (± 0.16) 


T p (MeV) 


149.59 ± 0.67 (± 0.02) 


(MeV/c 2 ) 


781.91 ± 0.18 (± 0.16) 


P^ (MeV) 


8.13 ± 0.36 (± 0.27) 


M 


(1.644 ± 0.061) x 10~ 3 


(rad) 


-0.011 ± 0.037 


m , (MeV/r 2 ) 


1493 + 15 


T p , (MeV) 


427 ± 31 




0.158 ± 0.018 


4> p , (rad) 


3.76 ± 0.10 


m p „ (MeV/c 2 ) 


1861 ± 17 


T pl , (MeV) 


316 ± 26 


Ivl 


0.068 ± 0.009 


<t>p« (rad) 


1.39 ± 0.20 


m p m (MeV/c 2 ) 


2254 ± 22 


I> (MeV) 


109 ± 76 


l<V"l 


n nn^i + 0034 

U.UUOI _ .0019 


<t> p w (rad) 


0.70 ± 0.51 



tween 0.9 and 1.2 GeV. This shows that the GS function, 
the parameters of which are mainly determined in the p 
peak region, together with the constraint at s' — 0, does 
not accurately describe the resonance tails. At higher 
masses, the validity of the VDM description that involves 
the parametrization of very broad resonances with large 
inelasticity, is somewhat arguable. However the overall 
agreement is satisfactory, notably in the 0.5-1.0 GeV re- 
gion. 

We compare the results for the resonance parameters 
(Table [Vl]) to those obtained by other experiments, not- 
ing that the comparison can be biased if the mass range 
or the parametrizations are different. The fitted p pa- 
rameters are compared to the results from CMD-2 [351 ] 
and SND [36]: for the mass m p , these two experiments 
obtain (776.0 ± 0.8) MeV/c 2 and (774.6 ± 0.6) MeV/c 2 , re- 
spectively, while for the width T p , they obtain (146.0 ± 
0.9) MeV and (146.1 ± 1.7) MeV. The fitted value of the 
phase 4>u] of the p — uj interference is not in good agree- 
ment with the CMD-2 value (0.182 ± 0.067 rad); SND 
uses a different parametrization. However in the CMD-2 
fit, the uj mass is fixed to the world-average value m u = 
782.65 MeV/c 2 [24j. If we fix m w to this value in the 
BABAR fit, the phase comes out to be (0.137 ± 0.023) rad, 
in agreement with CMD-2. In fact in the 18-parameter 
fit the fitted values for m w , and c u are strongly cor- 
related (80%). The fitted uj width r w is found to be con- 
sistent with the world-average value (8.49 ± 0.08) MeV 



obtained from the dominant 7r + 7r~7r° decay mode [24| . 

As the CMD-2 and SND experiments at Novosibirsk 
are well calibrated in energy with the resonant depolar- 
ization method, one can use the VDM fit to check the 
mass calibration by leaving the ui mass free, and using 
the CMD-2 result for the p — uj phase. One obtains 

m u = (782.68 ± 0.12 ± 0.27) MeV/c 2 , (34) 

where the first error is from the fit to the data and the 
second from the uncertainty on the CMD-2 value for </> w . 
The absolute difference with the world average uj mass is 

m u - m l UG = (0-03 ± 0.29) MeV/c 2 , (35) 

consistent with the calibration from the J/ifj study re- 
ported in SectionEm! (-0.16 ± 0.16) MeV/c 2 . 

J. Comparison to other experiments 

1. Pion form factor from e + e~ — > tt + it~ cross section 

The measured form factor FtAs') is compared to pub- 
lished data from the CMD-2 [3j|, SND and KLOE 
experiments. While the Novosibirsk results are obtained 
in the scan mode at fixed energy points, KLOE, like 
BABAR uses the ISR method, albeit at a much smaller 
energy (-y/s = 1.02 GeV). The KLOE data are ob- 
tained without direct detection of the ISR photon. More 
recently, KLOE has performed a new analysis |38[ where 
the ISR photon is detected at large angles, allowing them 
to collect data down to the threshold region. 

The data of the other experiments are compared with 
the result from the BABAR form factor fit, which was 
shown in the previous section to describe well the BABAR 
data itself. Each plot shows the relative difference be- 
tween the form factor-squared of the other experiment 
and BABAR as data points, while the width of the band 
around zero is the result of the propagation of statistical 
errors from the BABAR fit with systematic uncertainties 
in each mass region (Table [Vj added quadratically. 

The comparisons with other experiments are shown in 
Figs. [48] and [50l where the errors on the data points in- 
clude both statistical and systematic uncertainties. The 
agreement looks rather reasonable with the CMD-2 and 
SND measurements within systematic errors, the BABAR 
results lying generally above on the lower side of the p 
resonance. The discrepancy is larger with KLOE on and 
above the p peak. 

The region of the p — uj interference is examined in 
more detail in Fig. [49j No evidence is found for a signifi- 
cant variation in the steep part of the interference pattern 
around the uj mass, showing that the BAR4R mass cali- 
bration is not shifted with respect to CMD-2 and SND 
by more than 0.3MeV/c 2 . 

The comparison of the form factor-squared in the low- 
mass region is made in Fig. 1511 The agreement is reason- 
able, with some systematic departure with respect to the 
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FIG. 48: The relative difference of pion form factor-squared from the BABAR. fit in the 0.5-1 GeV region with CMD-2 (left) 
and with SND (right). Systematic and statistical uncertainties are included in the data points. The width of the BABAR band 
shows the propagation of statistical errors from the fit and the quoted systematic uncertainties, added quadratically. 
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FIG. 49: The relative difference of pion form factor-squared from the BABAR fit in the p — uj mass region with CMD-2 (left) 
and with SND (right). Systematic and statistical uncertainties are included in the data points. The width of the BABAR band 
shows the propagation of statistical errors from the fit and the quoted systematic uncertainties, added quadratically. 
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FIG. 50: The relative difference of pion form factor-squared 
from KLOE and the BABAR fit in the 0.5-1 GeV region. Sys- 
tematic and statistical uncertainties are included in the data 
points. The width of the BABAR band shows the propagation of 
statistical errors from the fit and the quoted systematic uncer- 
tainties, added quadratically. 
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FIG. 51: (color online). The measured pion form factor- 
squared compared to published results from other experiments. 
Systematic and statistical uncertainties are shown for all re- 
sults, with the diagonal elements of the BABAR, covariance ma- 
trix. 
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FIG. 52: (color online). The measured cross section for 
e + e~ — > 7r + 7r~(7) compared to published results from CMD-2 
up to 1.4 GeV and DM2 above. Systematic and statistical un- 
certainties are shown for all results, with the diagonal elements 
of the BABAR covariance matrix. 
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FIG. 53: The relative difference of pion form factor-squared 
from CMD-2 and SND and the BABAR fit in the region below 
0.5 GeV. Systematic and statistical uncertainties are included 
in the data points. The width of the BABAR band shows the 
propagation of statistical errors from the fit and the quoted 
systematic uncertainties, added quadratically. 
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FIG. 54: The relative difference of pion form factor-squared 
from CMD-2 and the BABAR fit in the region above 1 GeV. Sys- 
tematic and statistical uncertainties are included in the data 
points. The width of the BABAR band shows the propagation 
of statistical errors from the fit and the quoted systematic un- 
certainties, added quadratically. 
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NA7 experiment at CERN. A direct cross section com- 
parison is made in the large mass region in Fig. 1521 The 
BABAR results agree with CMD-2 up to 1.4 GeV, while 
the DM2 cross section [28[ appears to be larger by about 
30-40%. 

The comparison in relative terms of BABAR to other 
experiments is presented in Figs. [53] and [54] for masses 
lower than 0.5 GeV and between 1.0 and 1.4 GeV, respec- 
tively. The small discrepancy noticed between the BABAR 
fit and CMD-2 is in fact also observed in Fig. [47] where 
BABAR data are compared to the fit. So it points to a 
problem in the VDM parametrization rather than in the 
data. 



2. t spectral functions 



It is also appropriate to compare the present results 
to the r — » v T irir° spectral function. Taking isospin- 
breaking (IB) into account, the conserved vector current 



(CVC) relation between the e 
section with FSR a 



+, 



o 



" (7) bare cross 



(7) 

invariant mass distribution in r 
ified (ai |U as follows: 



and the normalized hadronic 
decays is mod- 
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where 
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27ra(0) 2 m2 
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R m (s) 



1 fp \ 3 \F (s)\i 
Gem(s) \P-J |F_(s)| 2 



(37) 



(38) 



Stttt and B e are the branching fractions for r decay into 
the v T irir° and v T eV e final states. Gem(s) is the long- 
distance QED radiative correction and Sew the short- 
distance electroweak radiative correction. F Q (s) and 
F_(s) are the electromagnetic and weak form factors, 
while /?o and /?_ are the pion velocities in the 7r + 7r~ and 
irir center-of-mass systems, respectively. 

Isospin-breaking (IB) corrections have been recently 
re-evaluated @ , and a new t analysis for the muon g — 2 
presented, taking advantage of the Belle data. It updates 
the previous analysis 2] . The Gem factor takes also into 
account p — u> interference, and the charged and neutral p 
mass difference, a charged and neutral p width difference 
from radiative decays, and the m^± — m^a mass difference 
in the form factor [4l|. Other recent approaches to IB 
breakin g h ave been considered, often based on specific 
models gl El| . 

Using the results from Ref. Q the corrected r and 
BABAR data can be compared directly. This is achieved in 
Figs. E5] for the ALEPH 0, CLEO [H, and Belle HH 
experiments, in a manner similar to the e + e~ compar- 
isons. Here there is another uncertainty resulting from 



the IB theoretical corrections, corresponding roughly to 
a scale uncertainty of 0.3%. For this comparison the 
spectral functions are normalized by the B xir = B(t — > 
TTTr^Vr) value measured by each experiment, rather than 
using the world average as usually done. In this way the 
spectral functions are really independent. The errors on 
the r data points include all sources of statistical and 
systematic uncertainties, B n7r and IB corrections. 

The comparison with ALEPH shows consistency 
within the systematic uncertainties up to the p peak and 
some slope above, with the remark that the ALEPH 
points are strongly correlated. Agreement is also ob- 
served within errors with the results of CLEO and Belle, 
CLEO being somewhat in between ALEPH and Belle 
above 0.8 GeV. 



X. THE 7T7T CONTRIBUTION TO THE 
ANOMALOUS MUON MAGNETIC MOMENT 

A. The BABAR Result 

The lowest-order loop contribution of the 7T7r(7) inter- 
mediate state to the muon magnetic anomaly is given 
by El 



r(7),LO 



OO 

^JdsK( S )a^ M (s), (39) 

1 mi. 



where K(s) is the QED kernel j48|, 
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FIG. 55: The relative difference of the form factor-squared from the r data of ALEPH (top), CLEO (middle) and Belle (bottom) 
with respect to the e + e~ — > ir + ir~ BABAR measurements in the 0.5-1 GeV region. Systematic and statistical uncertainties are 
included in the data points. The width of the BABAR band shows the propagation of statistical errors from the fit and the quoted 
systematic uncertainties, added quadratically. The r data are normalized to the value of B^n measured by each experiment 
independently. 



with x — (1 — f3fj)/{\ + /3 M ) and /3 M the muon velocity. 

The integration is carried out numerically over the 
measured cross section per mass bins. The statistical and 
systematic errors are computed using the corresponding 
covariance matrices described in Sections IIX El and IPC Fl 

Several tests are performed. 

• When the integral is performed with the origi- 
nal 50MeV bins of ISR luminosity the result is 
(514.40 ± 2.54 ± 3.11) x lfr 10 in the range 0.3- 



1.8 GeV, while the value (513.54 ± 2.22 ± 3.11) x 
10~ 10 is obtained with the chosen sliding-bin 
method. The difference is consistent with the sta- 
tistical fluctuations of the luminosity in the 50 McV 
bins (Fig. IM)) and the K(s) kernel weighting effect 
in the a M integral. 

In the 0.5-1.0 GeV range one compares the results 
obtained with the l p central' and the l p tails' condi- 
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tions. The main difference is the \ 2 selection, which 
affects the background level, the % 2 efficiency, and 
the mass resolution, hence the performance of the 
unfolding. For the 0.5-1.0 GeV range, the result 
of the integration with the 'central' conditions is 
445.94 x 10~ 10 in 2MeV bins, and 446.56 x 10~ 10 
with the 'tails' conditions in 10 MeV bins. Thus 
the effect of different resolution and efficiencies 
has little effect on the integral. The difference of 
0.62 x 10 -10 between the two analyses is consistent 
with their estimated non-common systematic errors 
and non-common statistical errors, which induce an 
uncertainty on the integral of 1.8 x 10 -10 . 

The evaluation of the integral in the threshold region 
was performed in previous estimates [2] using a polyno- 
mial expansion in s' for the pion form factor, incorpo- 
rating constraints on the normalization that -^(0) = 1 
and that the derivative of the form factor at s' = 
be given by the known quadratic charge radius of the 
pion. This procedure also compensated for the relatively 
poorer quality of data in this region. The BABAR continu- 
ous low-mass data permit a direct evaluation, consistent 
with the constrained method. The very small contribu- 
tion (0.55 ± 0.01) x 10~ 10 between the 2n threshold and 
0.3 GeV is evaluated using the extrapolation of the con- 
strained fit to the data between 0.3-0.5 GeV. 



TABLE VII: Evaluation of a^ (7) ' ° using the BABAR data 
in different mass regions (see text for details). The first error 
is statistical and the second systematic. 



m^Tt ran§ 


;e ( GeV) 


C h) ' LO (xlO- 10 ) 


0.28- 


-0.30 


0.55 ±0.01 ±0.01 


0.30- 


-0.50 


57.62 ±0.63 ±0.55 


0.50- 


-1.00 


445.94 ±2.10 ±2.51 


1.00- 


-1.80 


9.97 ±0.10 ±0.09 


0.28- 


-1.80 


514.09 ±2.22 ±3.11 



The BABAR results are given in Table IVIII in different 
mass ranges. The upper boundary (1.8 GeV) is chosen 
in accordance with previous evaluations [2j , in which the 
contribution of the higher energy region was computed 
using QCD. The contribution in the 1.8-3 GeV range, 
obtained with the present BABAR data, is indeed only 
(0.21 ± 0.01) x 10~ 10 , thus negligible with respect to the 
uncertainty in the main region. The contribution from 
threshold to 1.8 GeV is obtained for the first time from a 
single experiment: 

a ™-( 7 ),LO = (514.09 ± 2.22 ± 3.11) x 10" 10 , (41) 

where the first error is statistical and the second system- 
atic. The total uncertainty is 3.82 x 10 -10 , so that the 
precision of the measurement is 0.74%. 



B. Comparison to other determinations using e e 
and r data 

Direct comparison to the results from other experi- 
ments is complicated by two facts: (i) e + e~ scan experi- 
ments provide cross section measurements at discrete and 
unequally spaced energy values, while the ISR method 
provides a continuous spectrum, (ii) unlike BABAR, other 
experiments do not all cover the complete mass spectrum 
from threshold up to energies where the contributions 
become negligible. The latter problem is alleviated by 
appropriately combining different sets of measurements 
performed by the same experiment. Where gaps remain, 
they are filled by using the weighted-average cross sec- 
tion values from the other experiments. This approach 
has been thoroughly treated in Ref. [6], from which we 
extract the relevant integrals. The fraction of contribu- 
tions to the integrals estimated in this way ranges from 
3% for KLOE to 11% for CMD-2 (7% for SND), engen- 
dering some correlations between the total values (given 
in the energy range from 2m„ to 1.8 GeV). 

Correlations between systematic uncertainties have 
also been taken into account in Ref. particularly for 
radiative corrections, when combining the results from all 
experiments. The combination is achieved locally at the 
cross section level, taking into account possible disagree- 
ments leading to an increased uncertainty of the resulting 
average. The results are summarized in Table I Villi and 
Fig. [SB] and allow a direct comparison of all the deter- 
minations. They are indeed consistent within the errors, 
BABAR and CMD-2 being almost a factor of two more 
precise than SND and KLOE. Internal discrepancies be- 
tween the measurements increase the final uncertainty to 
3.2 x 10~ 10 , whereas the ideal value would be 2.4 x 10~ 10 
for fully consistent experiments. 

Similarly, the EABAR result is compared to determina- 
tions using r decays with IB corrections @ in Tabic IVTlTl 
and Fig. [S5J The agreement is found to be satisfactory, 
so the BABAR data reduces the previous tension between 
the e + e~ and r's values @. Regarding the consistency 
of all results, it is notable that the four inputs (CMD- 
2/SND, KLOE, BABAR, t) have completely independent 
systematic uncertainties. 



C. Impact of this result on the comparison of the 
Standard Model prediction and the direct 
measurement of the muon magnetic anomaly 

Even though the 2ir contribution is the dominant part 
of the hadronic LO VP component in the Standard Model 
prediction for the muon magnetic anomaly, all the con- 
tributions must be evaluated. The BABAR experiment 
has measured most of the relevant cross sections by the 
ISR method. Except for a few channels still unmeasured, 
BABAR results bring a new level of precision and domi- 
nate the picture. We follow here the recent analysis of 
Ref. I49| that uses all these measurements and those from 
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TABLE VIII: Evaluation of LO hadronic VP 2?r contribu- 
tions to the muon magnetic anomaly in the energy range 
[2m,r, 1.8 GeV] from BABAR, other e + e~ experiments 0], and 
r experiments [e| (see text for details). The errors are from 
both statistical and systematic sources. For the values derived 
from t decays, a common systematic error of 1.9 is included 
to account for uncertainties in the isospin-breaking correc- 
tions. Note that the combined results are not the weighted 
average of the different values, but originate from a proper 
local combination of the respective spectral functions @, |g. 



Experiment 


a ^( 7 ),LO (xl0 _l 0) 


BABAR 


514.1 ± 3.8 


CMD-2 


506.6 ±3.9 


SND 


505.1 ±6.7 


KLOE 


503.1 ±7.1 


combined e e~ 


507.8 ±3.2 


ALEPH 


508.7 ± 5.9 ± 1.9 


CLEO 


514.2 ± 10.4 ± 1.9 


OPAL 


526.9 ± 12.3 ± 1.9 


Belle 


513.7 ± 8.2 ± 1.9 


combined r 


515.2 ± 3.0 ± 1.9 
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FIG. 56: The LO hadronic VP 2n contributions to the muon 
magnetic anomaly, evaluated in the 2m 7r — 1.8 GeV range from 
the present analysis and other analyses using e + e~ data @] 
and t data [jj. The vertical bands show the values obtained 
for each data set by combining the respective spectral func- 
tions. 



other experiments. More recently, similar results have 
been obtained |50T |. 

Adding all contributions (QED, electroweak, hadronic 
LO VP other than 2ir, hadronic higher-order VP, 
hadronic light- by-light [51]) to the present 2n result 
alone, one obtains the predicted value of the muon mag- 



netic anomaly: 

a M = (11 659 186.5 ±5.4) x 10" 10 , (42) 

to be compared to the direct measurement [l[, slightly 
updated [521 ]: 

a° xp = (11 659 208.9 ±6.3) x 10" 10 . (43) 

The experimental value exceeds the prediction by 
(22.4±8.3) x 10~ 10 , i.e., 2.7 standard deviations. When 
the present BABAR cross section is combined with all 
available 2-7T data from e + e~ experiments the de- 
viation increases to (28.7 ± 8.0) x 10~ 10 , i.e., 3.6 stan- 
dard deviations. Although the deviation is not significant 
enough to claim a departure from the Standard Model, it 
confirms the trend of earlier results using previous e + e~ 
data 0-i] . 

XI. CONCLUSION 

The cross sections for the processes e + e~ — > 
M + / i_ (7)7iSR and e + e~ — » 7r + 7r~ (7)7131*, have been mea- 
sured by the BABAR experiment, where the additional 
photon may be produced either by FSR or ISR. Thanks 
to the properties of the ISR method, the corresponding 
e + e~ — > /i + /i _ (7) and e + e~ — > 7r + 7r~(7) cross sections 
have been determined from their thresholds to 3 GeV, 
thus covering completely the interesting region for calcu- 
lating hadronic vacuum polarization in the 7r + 7r~ chan- 
nel. 

For e + e _ — > p + /i~("f) the cross section is measured 
using the e + e~ luminosity, while the cross section for 
e + e~ — > 7r + 7r~(7), for which the highest precision is 
required, is obtained from the ratio of the radiative 
7r + 7r~ (7)7isr to /i + M~(7)7iSR mass spectra. In this way 
the pion results are independent of the e + e~ luminos- 
ity and important systematic effects cancel. As a major 
asset of the method, the pion-pair cross section is not 
sensitive to the model of radiative corrections in the gen- 
erator used for MC simulation. 

The measured absolute e + e~ — > (7)tisr cross 

section is found to agree with QED at NLO from thresh- 
old to 3 GeV, with a precision of 1.1% dominated by the 
e + e~ luminosity determination. 

The cross section for e + e~ — > 7r + 7r~(7) is obtained 
for the first time continuously in the energy range from 
threshold to 3 GeV. Its precision exceeds that of previous 
experiments in most of this range. The achieved system- 
atic uncertainty is 0.5% in the dominant p region from 
0.6 to 0.9 GeV. 

Fits of the pion form factor have been performed using 
a sum of contributions from isovector vector mesons: be- 
sides the dominant p resonance and isospin- violating p—uj 
interference, three higher states are needed to reproduce 
the structures observed in the measured spectrum. 

The results are in fair agreement with previous data 
from CMD-2 and SND, but some discrepancies are ob- 
served when compared to the KLOE data, particularly 
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on the p peak (3%) and above (up to 6% at 0.95 GeV). 
These differences exceed the uncertainties quoted by ei- 
ther experiment. The BABAR results are in agreement 
with the spectral functions derived from r — > ^ T 7r _ 7r° 
data, although some local deviations are seen in the line- 
shape at the 2% level with Belle. 

Finally, the BABAR results are used as input to the dis- 
persion integral yielding the 7r + 7r~(7) vacuum polariza- 
tion contribution at LO to the muon magnetic anomaly. 
This contribution amounts to (514.1±2.2 s t ot ±3.1 SJ/st ) x 
10~ 10 , the most precise result yet from a single experi- 
ment. This result brings the contribution estimated from 
all e + e~ — > Tt + n~("/) data combined in better agreement 
with the r estimate. When adding all other Standard 
Model contributions to the present 2tt result, in particu- 
lar using all available BABAR data on multihadronic pro- 
cesses, the predicted muon magnetic anomaly is found to 
be (11 659 186.5 ±5.4) x 1CT 10 , which is smaller than the 
direct measurement at BNL by 2.7<r. Adding all previous 
2ir data increases the deviation to 3.6cr. 

A claim for a breakdown of the Standard Model 
requires even more precise data, both for the direct 
anomaly measurement and the hadronic cross sections. 
But since a deviation of the observed size could be me- 
diated by new physics at a scale of a few 100 GeV, the 
present direct exploration for new phenomena performed 



at the Tevatron and the LHC will certainly bring valuable 
and complementary information. 
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